U  /D-A184  986 

UNCLASSIFIE 

ADVANCED 

RESEARCH 

NIPR-86-J 

CONCEPTS  THEORV  ANNUAL 
LA8  WASHINGTON  DC  12 

70 

REPORT  1986<U)  NAVAL 

UG  87  NRL -HR-8851 

F/G  28/9 

NL 

■ 

■ 

r 

a 

■ 

a 

in 

y 

1 

■■ 

Kpl 

*w 

i 


a 


v.J 

y>, 

v • 


I 

ny 

i*  ^ 


tV:iJ 

iW 

% 

*# 


128  12.5 


125  U  H.6 


MICROCOPY  RESOLUTION  TEST  CHART 
NATIONAL  BUREAU  OF  STANDARDS  196.5  A 


M.  a a&  avail- n  ■*■■■ 


AD- A 184  986 


NRL  Memorandum  Report  6051 


Advanced  Concepts  Theory  Annual  Report  1986 

Plasma  Radiation  Branch 
Plasma  Physics  Division 


August  12,  1987 


This  research  was  sponsored  by  the  Defense  Nuclear  Agency  under  Subtask  RL  RA,  work  unit  00049, 
and  work  unit  title  “X-Ray  Source  Development  Theory,"  MIPR  No.  86-570. 


Approved  for  public  release;  distribution  unlimited. 

87  9  18  004 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 

2a.  SECURITY  CLASSIFICATION  AUTHORITY 


REPORT  DOCUMENTATION  PAGE 

I'b  RESTRICTIVE  MARKINGS 


2b  DECLASSIFICATION  /  DOWNGRADING  SCHEDULE 
4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

NRL  Memorandum  Report  6051 

6a.  NAME  OF  PERFORMING  ORGANIZATION  (6b 

Naval  Research  Laboratory  _ 

6c  ADDRESS  (City.  State,  and  ZIP  Code) 

Washington.  DC  20375-5000 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution  unlimited. 

5  MONITORING  ORGANIZATION  REPORT  NUM8£R(S) 


6b  OFFICE  SYMBOL  7a  NAME  OF  MONITORING  ORGANIZATION 
(If  applicable) 

Code  4720  Defense  Nuclear  Agency 


7b  ADDRESS  (City,  State,  and  ZIP  Code) 

6801  Telegraph  Road 
Alexandria,  VA  22310 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Defense  Nuclear  Agency 

8c.  ADDRESS  (City,  state,  and  ZIP  Code) 

Alexandria.  VA  22310 


1 1  TITLE  (Include  Security  Classification) 

Advanced  Concepts  Theory  Annual  Report  1986 

12  PERSONAL  AUTHOR(S) 

Plasma  Radiation  Branch 


8b.  OFFICE  SYMBOL  I  9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 
(If  applicable)  j 


10  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM  PROJECT 

ELEMENT  NO  NO 

627 15H 


PROJECT 

TASK 

NO 

NO. 

WORK  UNIT 
ACCESSION  NO 

DN880-191 


14  DATE  OF  REPORT  (Year,  Month.  Day)  h  5  PAGE  COUNT 

12  August  1987  I  108 


13a.  TYPE  OF  REPORT  1 3b  TIME  COVERED  14  DATE  OF  REPORT  (Year,  Month.  Day)  15  PAGE  COUNT 

Annual  from  1/86  to  12/86  12  August  1987  _ 108  _ 

16  SUPPLEMENTARY  NOTATION  .  ,  .  ,  r  . 

This  research  was  sponsored  by  the  Defense  Nuclear  Agency 
under  Subtask  RL  RA.  work  unit  00049,  and  work  unit  title  “X-Ray  Source  Development  Theory."  MIPR  No.  86-570. 

17  _ _ COSATI  COOES _  18  SUBJECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

field  group  sub-group  Z  pinch  physics  Radiation  hydrodynamics 

Plasma  mixtures  Gyrokinetic  propagation 

Plasma  opening  switches _ 

19  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number) 

A  This  report  details  the  work  of  the  Plasma  Radiation  Branch  conducted  in  FY  86  with  respect  to  several 
critical  problems  concerning  the  design  of  Z  pinch  implosions  and  the  optimization  of  their  radiative  output. 

Separate  sections  describe  progress  in  atomic  modeling  of  K  shell  copper  radiation,  radiation  hydrodynamics, 
argon  implosions,  design  of  anti-zippering  nozzles,  and  analyses  of  anomalous  resistivity  and  plasma  opening 
switches. 

'V 


17 

COSATI  CODES  1 

FIELD 

GROUP 

SUB-GROUP 

20  DISTRIBUTION/ AVAILABILITY  OF  ABSTRACT  21  ABSTRACT  SECURITY  CLASSIFICATION 

E3  UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT  □  OTIC  USERS  UNCLASSIFIED _ 

22a  NAME  OF  RESPONSIBLE  INDIVIDUAL  22b  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 

Jack  Davis  (202)  767-3278  Code  4720 

DD  FORM  1 473,  84  MAR  83  APR  edition  may  be  used  until  exhausted  SECURITY  CLASSIFICATION  OF  THIS  PAGE 

All  other  editions  are  obsolete 


CONTENTS 


EXECUTIVE  SUMMARY  . 

A.  COPPER  ATOMIC  MODEL  AND  APPLICATION  TO  HARD 

Z-PINCH  X-RAYS  . 

B.  KRYPTON-TUNGSTEN  MIXTURES  IN  Z-PINCH  PLASMA  IMPLOSIONS  .. 

C.  K-SHELL  YIELD  FROM  AN  ARGON  GAS  PUFF  . 

D.  ANTI-ZIPPERING  NOZZLE  IMPLOSION  CHARACTERISTICS  AND  THEIR 

IMPORTANCE  TO  UNIFORMITY  IN  A  Z-PINCH  . 

E.  OPENING  SWITCH  AND  DIODE  MODELING  . 

F.  APPENDIX  —  EFFECT  OF  PULSE  SHARPENING  ON  IMPLODING 

NEON  Z-PINCH  PLASMAS  . 


Accession  For 

NTIS  GRA&I 

J? 

DTIC  TAB 

Unannounced 

□ 

Justification _ 

By - 

Distribution/ _ 

A  _ _ 


HJ 


11 

pe 


[»I«I 


Executive  Summary 


During  FY  86  the  Plasma  Radiation  Branch  at  NRL  has  focused  its 
resources  on  several  critical  problems  concerning  the  design  of  Z  pinches 
and  the  fundamental  physical  processes  governing  the  evolution  of  their 
implosions  and  radiative  output.  Our  1986  Advanced  Concepts  Theory  Annual 
Report  consists  of  five  self-contained  sections,  each  reporting  substantial 
progress  on  significant  design  innovations  or  Z  pinch  physics. 

The  first  two  of  these  sections  report  progress  on  atomic/hydrodynamic 
modeling  of  copper  and  krypton  implosions.  Achieving  substantial  K  shell 
emission  with  either  or  both  of  these  high  atomic  number  elements  would 
represent  a  quantum  leap  in  DNA’s  simulation  capability.  The  K  shell  lines 
in  these  elements  range  in  energy  from  8  to  15  keV.  Section  A  reports  on 
the  development  of  a  detailed  K  shell  atomic  model  for  copper  ions  and  the 
benchmarking  of  the  spectra  generated  by  this  theoretical  model  against 
published  experimental  spectra.  The  temperature-density  diagnostics 
available  from  the  copper  K  shell  spectrum  have  been  identified  and  in 
FY  87  will  be  applied  to  the  copper  wire  experiments  conducted  on  BLACKJACK 
5  at  Maxwell  Laboratories,  Inc. 

Section  B  explores  a  concept  which  has  appeared  frequently  in  the  open 
Soviet  literature  and  is  considered  by  many  Soviet  scientists  to  be  an 
established  fact  -  enhanced  compressions  of  high  atomic  number  plasmas 
driven  by  radiative  cooling.  The  idea  is  that  very  high  densities  and 
hence  radiative  yields  are  obtainable  by  seeding  a  principal  radiating 
element  (krypton  gas  puff,  in  our  case)  with  a  very  effective  high  atomic 
number  "coolant"  (we  have  chosen  tungsten).  As  the  tungsten  coolant 
radiates,  the  temperature  and  pressure  of  the  plasma  initially  remain  lower 
than  would  be  the  case  with  a  pure  krypton  plasma.  This  permits  a  greater 
degree  of  compression  to  be  achieved,  and,  ultimately,  increases  the 
radiative  yield  which  scales  as  the  square  of  the  density.  As  the  results 
of  Section  B  demonstrate,  the  concept  appears  to  be  promising  and  we 
strongly  recommend  that  DNA  invest  in  additional  modeling  (in  particular 
improvement  of  our  high  Z  atomic  data  base)  and  some  preliminary 
experiments  both  with  and  without  high  Z  coolants  included  to  provide  a 
laboratory  test  of  this  yield-enhancing  design. 


Section  C  presents  detailed  one-dimensional  self-consistent  radiative- 
hydrodynamic  modeling  of  argon  implosions  on  Gamble-II.  It  is  well  known 
that  as  the  mass  injected  into  a  gas  puff  pinch  is  increased,  the  radiative 
yield  at  first  increases,  but  reaches  a  maximum  and  then  falls 
substantially  at  high  mass  loads.  The  analysis  of  this  section  identifies 
and  explains  the  physics  traoeoffs  which  result  in  this  phenomenon.  Even 
though  more  kinetic  energy  can  be  stored  in  a  lighter  mass,  there  may  not 
be  enough  radiating  ions  present  to  effectively  convert  this  energy. 
Conversely,  at  high  mass,  the  smaller  kinetic  energy  which  can  be  imparted 
to  a  heavy  load  prevents  high  yields  from  being  achieved.  When  the  fields 
are  coupled  to  the  plasma  near  the  time  of  peak  current,  full  two- 
dimensional  code  analysis  is  required  to  simulate  the  formation  of  the 
various  "hot  spots"  arising  from  plasma  instabilities.  Nevertheless,  the 
1-D  simulations  have  successfully  predicted  within  experimental  limits  the 
mass  loading  point  at  which  Gamble-II '3  radiation  yields  maximize. 

One  of  the  most  important  deleterious  effects  limiting  the  peak  power 
achievable  by  Z  pinch  implosions  is  known  as  "zippering"  -  the  non  uniform 
pinch  assembly  along  the  Z  axis  which  spreads  the  radiation  pulse  in 
time.  Frequently  the  total  K-shell  pulse  is  nearly  an  order  of  magnitude 
longer  than  would  be  the  case  if  the  entire  length  of  plasma  assembled 
simultaneously.  Section  D  presents  results  from  hydrodynamic  modeling  of 
various  nozzle  designs  intended  to  counteract  the  zippering  effect  by 
tilting  the  nozzle  to  spray  the  cold  gas  inward.  Results  show  that  certain 
nozzle  designs  can  result  in  the  simultaneous  assembly  of  a  length  of  2  cm 
in  the  Z  direction.  The  next  step  in  this  investigation  is  the  direct 
simulation  of  radiative  power  pulse  shapes,  using  anti-zipper ing  nozzles. 

The  scaling  of  radiative  yields  has  been  observed  to  follow 
approximately  the  fourth  power  of  peak  current  (see,  for  instance,  the 
paper  in  the  Appendix).  This  implies  that  a  detailed  knowledge  of  plasma 
resistivity  and  its  effect  on  current  distribution  could  yield  impressive 
dividends  in  design  of  machines  and  loads  to  optimize  radiation.  For 
instance,  a  20%  increase  in  current  would  boost  the  yield  by  a  factor  of 
two.  Unfortunately ,  the  theory  of  anomalous  resistivity  in  plasmas  is  not 
yet  firmly  established  for  Z  pinch  geometries.  However,  our  gyrokinetic 
propagation  studies  (Section  E)  are  providing  an  increasingly  firm  base  for 
anomalous  resistivity  calculations  in  radiation  source  plasmas  (PRS's)  and 


in  opening  switch  plasmas  (VOS' s).  In  VOS  schemes,  the  power  flow  in  the 
tenuous  magnetized  regions  will  influence  the  energy  transfer  efficiency 
opening  times  of  the  devices. 

Simulations  of  plasma  opening  switches  are  confirming  the  low  density, 
not  electron  mode  of  operation.  Mechanisms  that  produce  a  wider  current 
channel  have  been  identified.  Simple  electrostatic  diode  simulations  have 
proven  adequate  for  treatment  of  realistic  geometries  such  as  ring  diodes 
and  a  coordinated  effort  to  address  design  details  is  proceeding  with 
Maxwell  Laboratories. 

Finally,  the  Appendix  contains  the  final  version  of  a  paper  "Effect  of 
Pulse  Sharpening  on  Imploding  Neon  Z^Pinch  Plasmas"  which  describes  the 
greater  radiating  efficiency  of  Gamble^II  when  equipped  with  a  plasma 
erosion  opening  switch.  Some  of  this  material  appeared  in  different  form 
in  last  year's  Advanced  Concepts  Final  Report  (NRL  Memorandum  Report 
5759).  The  published  paper  is  included  here  for  convenient  reference. 


In  addition  to  the  work  presented  here  there  have  been  a  number  of 
reports  on  related  activities.  They  include: 

1 .  NRL  Memorandum  Report  5765  "Gamble-II  Imploding  Sodium  Plasma.  I 
Calibration  of  the  Heliumlike  Resonance  Line  as  a  Pump  Source  and 
Detection  of  Fluorescence  in  Neon”,  J.  Davis,  J.  E.  Rogerson,  and  J.  P. 
Apruzese  ,  10  April  1986. 

2.  NRL  Memorandum  Report  5771  "Argon  Puff  Gas  Soft  X-Ray  Laser",  J. 

Davis,  J.  P.  Apruzese,  C.  Agritellis  and  P.  Kepple,  April  1986. 

3.  NRL  Memorandum  Report  5776  "Gamble-II  Imploding  Sodium  Plasma.  II. 
Uniformly  Filled  Z-Pinch",  J.  Davis,  J.  E.  Rogerson,  and  J.  P. 

Apruzese,  May  1986. 

4.  NRL  Memorandum  Report  5871  "A  Comparison  of  Radiation  Transport  Methods 
in  Axisymmetric  Geometries",  J.  Giuliani,  October  1986. 

5.  "Dynamics  of  Imploding  Neon  Gas-Puff  Plasmas",  R.  W.  Clark,  J.  Davis 
and  F.  L.  Cochran,  Phys.  of  Fluids  Vol.  29  (6)  p.  1971,  June  1986. 

6.  "X-Ray  Lasing  in  a  Na/Ne  Plasma  Environment",  F.  L.  Cochran,  J.  Davis 
and  J.  P.  Apruzese,  Journal  of  Applied  Physics  5J_  p.  27,  January  1985. 
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ADVANCED  CONCEPTS  THEORY  ANNUAL  REPORT  1986 


Section  A.  Copper  Atomic  Model  and  Application  to  hard  Z- 
Pinch  X-Rays 


It  has  long  been  considered  highly  desirable  by  DNA  to 
obtain  a  substantial  yield  of  x-rays  in  the  10  keV  region. 
Copper's  chemical  and  mechanical  properties  facilitate 
target  fabrication  and  it  is,  therefore,  an  excellent  choice 
of  material  for  exploratory  experiments  aimed  at  producing 
copious  K-shell  radiation.  We  have  studied  the  ionization 
to  the  K-shell  in  copper  and  the  resulting  high  energy 
radiation  as  a  function  of  radius,  density  and  temperature 
for  cylindrical  geometries  using  a  collisional  radiative 
equilibrium  ( CRE )  model.  A  simple  K-shell  copper  model, 
which  has  been  benchmarked  against  experiment1,  is  used  to 
analyze  the  feasibility  of  employing  copper  in  Z-pinches  to 
produce  significant  amounts  of  K-shell  radiation. 


I.  MODEL  DESCRIPTION 


Since  we  consider  copper  plasmas  of  only  such  temperatures 
and  densities  that  the  K-shell  dominates  the  ionic  species, 
we  employ  a  multistate,  multilevel  model  which  includes 


1  The  krypton  rates  were  downscaled  to  copper  using  a 
computer  code,  which  has  since  been  modified  so  that  the 
atomic  rates  from  any  element  can  be  hydrogeni cal ly 
downscaled  to  those  of  any  other  element.  This  means  that 
rate  tables  (covering  the  K-shell  region)  can  be 
constructed  fairly  quickly  for  different  elements. 

This  code  will  be  a  valuable  tool  for  DNA  in  analyzing 
the  feasibility  of  different  elements,  or  combination  of 
elements,  as  plasma  sources  for  Z-pinches. 
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excited  levels  through  principal  quantum  number  n-5  for 
lithiumlike,  heliumlike,  and  hydrogenlike  copper  in  addition 
to  all  thirty  ground  states  (neutral  through  fully 
stripped).  For  lithiumlike  Cu  XXVII,  the  six  excited  states 
included  are:  ls^2p  ^P,  ls^3s  ^S,  ls^3p  ^P,  ls^3d  ,  n=4 
doublet,  and  n=5  doublet.  For  heliumlike  Cu  XXVIII,  there 
are  eight  excited  states:  ls2s  JS,  ls2s  AS,  ls2p  JP, 
ls2p  ^P,  the  consolidated  n=3  triplet  states,  the 
consolidated  n=3  singlet  states,  and  the  consolidated  n-4 
and  n=5  levels.  Finally,  for  hydrogenlike  Cu  XXIX  there  are 
four  excited  states:  the  n=2,  3,  4  and  5  levels,  with  the 
sublevels  treated  as  statistically  populated. 

The  copper  rates  were  obtained  by  hydrogenically 
downscaling  krypton  transition  rates. ^  The  krypton  rates 
themselves,  however,  were  not  calculated  hydrogenically. 

For  details  of  the  krypton  model,  the  reader  may  consult 
Ref.  2  which  describes  an  argon  model  similar  to  the  krypton 
model  used  here.  The  atomic  processes  populating  and 
depopulating  the  levels  are:  collisional  ionization, 
collisional  and  radiative  recombination,  collisional 
excitation  and  deexcitation,  dielectronic  recombination  and 
spontaneous  radiative  decay. 

The  rate  equations  are  solved  in  the  steady-state 
approximation  for  a  homogeneous  plasma  of  constant  ion 
density  and  uniform  temperature.  An  iterative  technique^  is 
used  to  obtain  self-consistency  between  the  fractional  level 
populations  and  the  radiation  field.  To  account  for  opacity 
effects,  a  detailed  frequency  grid  spanning  the  region  of 
interest  (8  to  11  keV)  has  been  formed  by  combining  fine 
grids  centered  on  each  line  in  the  region.  From  8  to  11 
keV,  radiation  is  then  transported  at  each  grid  point; 
outside  of  this  region,  individual  lines  are  transported 

4 

probabilistically  . 

The  cylindrical  geometry  considered  here  is  a  1  mm 
diameter,  3  cm  long  cylindrical  plasma  whose  temperature  and 
ion  density  are  varied.  (Note  that  the  total  number  of  ions 
is  not  constant  in  this  case.)  This  geometry  constitutes  an 
approximation  of  Z-pinch  plasmas.  In  the  cylindrical  runs, 
the  ion  and  electron  densities  are  low  enough  that  line 
shapes  can  be  described  adequately  by  Voigt  profiles. 


II.  COMPARISON  WITH  EXPERIMENT 


Yaakobi  and  Burek^  have  published  a  spectrum  obtained  by 

imploding  a  thin  copper  shell.  They  imploded  a  0.3  n m  thick 

copper  shell  of  diameter  80  */m,  obtaining  the  time- 

integrated  spectrum  shown  in  Fig.  1.  There  was  film 

saturation,  so  that  the  He  a  line  is  actually  more  intense 

than  it  appears  to  be.  We  modeled  the  plasma  conditions  by 

assuming  that  the  spectrum  was  obtained  primarily  at  peak 

compression  and  that,  at  that  time,  the  plasma  was 

1 4 

homogeneous,  containing  a  total  of  5.14x10  copper  ions. 

By  varying  the  temperature  and  ion  density  (equivalent  to 

sphere  diameter  at  peak  compression)  we  obtained  a  best  fit, 

Fig.  2,  for  a  temperature  of  1.3  keV  and  ion  density  of 
22  -3 

3.64x10  cm  (or  30  yt/m  diameter).  Since  doubly  excited 
levels  have  not  been  included  in  the  atomic  model,  not  all 
of  the  lines  that  appear  in  the  experimental  spectrum  show 
up  in  Fig.  2.  According  to  our  calculations,  at  4  keV,  or 
even  at  2  keV,  the  hydrogenlike  Ly  a  line  at  8.727  keV  is 
much  too  intense  to  agree  with  the  experimental  spectrum. 
Indeed,  at  4  keV  the  He  a  and  Ly  a  are  of  compariable 
intensity,  whereas  in  Fig.  1,  the  Ly  a.  line  has  barely 
registered. 

III.  RESULTS 

In  the  cylindrical  case,  spectra,  fractional  populations 

and  He  a  line  power  outputs  were  calculated  for  a  1  mm 

diameter,  3  cm  long  homogeneous  plasma  at  various 

temperatures  and  densities.  The  densities  assumed  in  this 

19  -  3 

geometry  are  around  2x10  cm  ,  considerably  less  than 
typically  reached  in  a  spherical  implosion.  As  a 
consequence  there  is  less  collsional  excitation.  This  means 
that  there  will  be  fewer  ions  in  the  excited  heliumlike 
levels,  and  since  these  excited  levels  can  be  ionized  to 
produce  the  hydrogenlike  ions^,  there  will  be  fewer 
hydrogenlike  ions.  Figures  3  through  6  are  theoretical 
spectra  for  temperatures  of  1.5,  2,  4,  and  6  keV, 
respectively. 
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Figure  7  displays  hydrogenlike  and  heliumlike  copper 

fractions  as  a  function  of  temperature  for  the  plasma 

cylinder  described  above  at  ion  densities  ranging  from 
19  19  -  3 

2x10  '  to  5x10  cm  .  An  important  facet  of  Fig.  7  is  that 
while  heliumlike  copper  is  produced  in  substantial  abundance 
(45%)  at  a  plasma  temperature  of  2  keV,  the  hydrogenlike 
copper  does  not  appear  at  even  the  10%  level  until  4.5  keV. 
Therefore,  there  will  be  a  regime  in  which  only  the 
heliumlike  lines  are  observed.  The  relative  ease  in 
producing  the  heliumlike  stage  relative  to  the  hydrogenlike 
stage  is  a  straight  forward  consequence  of  the  ionization 
potentials.  It  requires  2.6  keV  to  ionize  a  lithiumlike 
copper  ion  to  heliumlike;  but  to  ionize  a  heliumlike  copper 
ion  to  the  hydrogenlike  stage  requires  11  keV. 

This  phenomenon  is  reflected  in  line  power  outputs.  Even 
at  a  temperature  of  4  keV,  the  hydrogenlike  line  power  is 
about  an  order  of  magnitude  below  the  heliumlike  line 
output.  We  have  fit  the  K-shell  power  output  in  J/nsec  to 
the  following  scaling  law,  good  to  10%  for  temperatures 
between  1.5  and  3.5  keV,  and  good  to  28%  for  temperatures  up 
to  5  keV: 

P ( ni , T )  =  130  nA2  exp[ -12/T2 ]  J/nsec, 

19 

where  n.  is  the  total  copper  ion  density  in  units  of  1x10 

_  3  1 

cm  ,  T  is  the  temperature  in  keV,  and  we  have  assumed  a  3 

cm  long  copper  plasma  of  diameter  1  mm.  At  an  ion  density 
19  -  3 

of  2x10  cm  ,  assuming  a  15  ns  radiating  pulse,  a  copper 
K-shell  yield  of  1  kJ  would  be  attained  at  a  plasma 
temperature  of  2.5  keV. 

Figure  8  shows  the  total  power  output  due  to  hydrogenlike, 

heliumlike,  and  lithiumlike  ions  for  1  cm  long  cylinders 

containing  10  and  100  u g  of  copper  at  2  keV.  At  large  radii 

2 

the  power  output  follows  the  1/r  law.  As  the  radius 

decreases  and  the  density  increases,  systematic  variations 

in  the  ionization  stages  (here  between  lithiumlike  and 

berylliumlike  ions)  cause  the  power  output  to  fall  below  the 
2 

1/r  line.  A  more  complete  model,  containing  excited  levels 

2 

of  beryllium,  would  follow  the  1/r  line  for  smaller  radii 
more  closely. 


IV.  CONCLUSION 


A  K-shel I  copper  model  has  been  used  to  examine  the 
ionization  and  spectral  output  of  Z-pinches  in  the  &  .  . 

keV  region.  For  the  densities  considered  here,  .*  is  r  * 
until  about  7  keV  that  the  hydrogenlike  population  is  ha.f 
that  of  the  heliumlike  fractional  population.  This  means 
that  for  temperatures  less  than  7  keV,  there  will  be  - n . y 
very  weak  hydrogenlike  lines.  A  large  fraction  of  the  ::rs 
tend  to  remain  in  the  heliumlike  ground  state  even  at  high 
temperatures,  causing  the  He  a  power  output  at  line  center 
to  be  very  intense.  Yields  of  a  kilojoule  are  obtainable  a" 
temperatures  of  2.5  keV.  The  relative  intensiites  of  the 
hydrogenlike,  heliumlike,  and  any  inner  shell  lines  present 
will  provide  valuable  clues  to  the  condition^  achieved  in 
the  forthcoming  copper  wire  experiments  at  Maxwell 
Laboratories . 
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FIG.  1.  Experimental 
and  Burek  (Ref.  1)  by 
thick  copper  shell  of 
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FIG.  3.  Copper  spectrum  obtained  with  collisional- 
radiative-equilibrium  plasma  model  for  a  3  cm  long,  1mm 
diameter  cylinder  at  a  temperature  of  1.5  keV  and  ion 
density  of  2x10  cm  .  (Electron  density  is  4.99x10  u 
cm-3  )  . 
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FIG.  4.  Same  as  Fig.  3,  except  the  temperature  is  2  keV 

20  -3 

(Electron  density  is  5.27x10  cm  ). 
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Section  B.  Krypton-Tungsten  Mixtures  in  Z-Pinch  Plasma  Implosions 


Existing  pulse  power  generators  have  not  yet  produced  large  amounts  of 
line  x-rays  in  the  5-10  keV  range.  Recently,  efforts  have  been  undertaken 
to  explore  ways  to  increase  the  production  of  line  x-rays  in  PRS  using  a 
variety  of  approaches.  The  copper-wire  implosion  approach  is  discussed  in 
Section  A  above.  Such  approach  involves  seeding  a  single  wire  or  gas  puff 
plasma  with  a  high-Z  tracer  to  determine  whether  the  radiative  cooling  due 
to  the  tracer  is  sufficient  to  increase  the  plasma  compression  in  a  z-pinch 
implosion  resulting  in  higher  densities  and  greater  radiative  yields  of 
high  energy  photons. 

In  this  study  the  implosion  of  a  cylindrical  annular  Krypton  (Kr)  gas 

puff  plasma  seeded  with  a  small  amount  of  tungsten  (V)  tracer  is  studied. 

The  theory  describing  the  implosion  is  based  on  a  dynamic  z-pinch  model 

self-consistently  coupled  to  a  Collisional  Radiative  Equilibrium  (CRE) 

ionization  dynamics  and  radiative  transport  model^.  This  model  provides  a 

self-consistent  picture  of  the  temperature,  density,  size,  and  level 

populations  during  the  plasma  evolution.  Any  pathological  behavior,  such 

as  instabilities,  is  ignored  in  this  simplified  approach;  however,  general 

trends  in  the  plasma  behavior  should  be  indicated  by  this  model,  and  it 

should  serve  as  a  useful  guide  for  more  detailed  research. 

The  principal  plasma  species  Kr  is  comprehensively  treated.  A  low 

2 

density  representation  is  used  for  the  V  impurity  .  Opacity  effects  for 

the  V  radiation  are  taken  into  account  in  an  ad  hoc  manner  by  calculating 

2 

the  ratio  of  total  emitted  radiation  to  coronal  radiation  for  Kr  ,  and 

scaling  down  the  U  coronal  radiation  accordingly. 

The  implosions  modeled  in  this  study  are  driven  by  high  current 

discharges,  which  are  based  on  data  obtained  from  Physics  International 

3 

Company  regarding  the  output  of  their  Double-EAGLE  generator  . 

The  first  current  profile  considered  was  the  triangular  shaped  pulse 
shown  in  Figure  1.  Comparisons  between  the  responses  of  a  pure  Kr  plasma 
with  a  mixture  of  90%  Kr  -  10%  V  by  mass  were  made.  Results  for  different 
initial  cylindrical  widths  for  the  plasma  were  calculated  for  mass  loadings 
of  175,  200,  and  225  ugm/cm.  The  cylinder  length  was  taken  to  be  3  cm, 
and  the  initial  outer  and  inner  radii  were  chosen  so  that  the  average 
plasma  radius  was  initially  about  1.5  cm.  Figures  2a  -  2d  shows  the 

variations  of  the  minimum  plasma  radius  and  the  maximum  ion  density, 

temperature,  and  total  radiation  with  initial  plasma  shell  width  and  mass. 


In  each  case,  the  V-seeded  plasma  gives  a  better  implosion  than  the  pure  Kr 
case.  In  this  model,  this  effect  is  due  to  the  enhanced  radiative  cooling 
from  the  W  impurity. 

A  second  set  of  calculations  was  performed  using  the  current  profile 
shown  in  Figure  3.  The  use  of  the  second  profile  allows  for  comparing  the 
effect  of  a  flat  peak  with  25  nsec  duration  with  the  initial  triangular 
shaped  profile.  For  these  simulations,  the  mass  was  kept  constant  at  225 
Wgm/cm,  but  calculations  for  pure  Kr,  90%  Kr  -  10%  W,  and  80%  Kr  -  20%  V 
were  performed;  Figure  4  summarizes  these  results  for  the  same  extremal 
parameters  shown  in  Figure  2.  These  curves  again  indicate  that  seeding  the 
Kr  plasma  with  a  V  impurity  enhances  the  implosion,  and  that  the  larger  the 
amount  of  the  impurity,  the  better  the  implosion. 

Comparison  of  Figures  2  and  4  for  the  100%  Kr  and  90%  Kr  -  10%  V 
results  indicate,  as  might  be  expected,  that  the  flat-topped  current  trace 
in  Figure  3  gives  lower  minimum  radius  and  higher  maximum  ion  densities  and 
total  radiation  than  the  current  pulse  in  Figure  1.  For  the  temperature, 
however,  this  trend  is  only  observed  at  the  higher  plasma  widths.  The 
reasons  for  this  temperature  behavior  are  not  clear.  Another  interesting 
feature  of  Figure  2  is  the  apparent  virtual  independence  of  maximum 
temperature  on  plasma  width,  whereas  in  Figure  4  there  is  a  stronger 
dependence. 

It  should  be  emphasized  that  the  numerical  results  indicated  on  these 
figures  are  obtained  from  a  zero-dimensional  simplified  treatment  of  the 
dynamics  and  should  be  interpreted  accordingly.  A  more  sophisticated 
treatment  would  involve  multi-zone  calculations  in  one  or  more  dimensions 
with  coupled  hydrodynamics  and  radiation  transport.  Also,  the  impurity 
element  should  be  treated  on  an  equal  footing  with  the  principal  species; 
i.e.,  using  a  full  atomic  model  with  radiation  transport  for  both.  Such  a 
program  could  become  very  costly  in  terms  of  computer  time  and  storage. 

CONCLUSIONS 

Seeding  with  trace  amounts  of  high-Z  elements  produces  the  desirable 
effects;  better  compressions,  higher  densities,  and  greater  radiative 
yields.  This  method  of  enhancing  yields  appears  very  promising  and  will  be 
investigated  further.  If  trace  amounts  of  one  element  is  promising,  trace 
amounts  of  several  elements  may  be  of  even  greater  utility  as  opacity 

4 

effects  on  the  radiative  cooling  may  be  reduced  in  this  faction. 
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Section  C.  K-Shel I  Yield  from  an  Argon  Gas  Puff 

This  section  explores  theoretically  the  determination  of  the  parameter 
space  under  which  Gamble  II  -  like  devices  produce  maximum  K-shel I 
emission.  In  particular  we  try  to  optimize  the  mass  loading  and  initial 
radial  position  of  the  gas  puff. 

The  work  is  done  by  computer  simulation  using  a  one-dimensional  non- 
LTE  radiation  hydrodynamics  model.  Because  the  impedance  of  the  plasma  can 
at  times  influence  the  electrical  characteristics  of  the  machine,  we  have 
decided  to  drive  the  system  with  an  open  circuit  voltage  profile  rather 
than  a  current  trace. 

a .  Mode  I 

The  above  parameter  study  is  carried  out  using  a  1-dimensional,  30 
zone,  2  temperature,  lagrangian  hydrodynamics  code  in  which  the  ionization 
and  radiation  dynamics  is  treated  self-consistently  using  the  collisional 
radiative  equilibrium  (CRE)  radiation  and  atomic  physics  models  developed 
and  used  by  the  NRL  Plasma  Radiation  Branch*  over  the  last  several  years. 

In  particular  the  argon  model  employed  utilizes  19  ground  states  ,  44 
excited  levels  and  56  bound-bound  transitions.  The  excited  states  were 
chosen  to  predominantly  surround  hydrogen- I i ke  through  lithium-like  argon 
in  order  to  provide  good  information  on  K-shel I  radiative  yields,  however 
excited  state  structure  is  included  for  ion  levels  up  to  sodium- 1  ike  argon. 

In  the  subsequent  studies  we  are  imploding  an  annulus  of  argon  gas 
that  is  .3  cm  between  inner  and  outer  radius.  The  radius  R  is  ini  tally 
fixed  at  1.35  cm  and  the  mass  density  is  varied  from  5  to  60  p gm/cm.  Once 
an  optimum  mass  loading  is  found  we  fix  the  mass  loading  at  this  value  and 
then  vary  the  initial  radial  position  from  .25  to  1.35  cm.  The  puff  is 
imploded  assuming  the  open  circuit  voltage  profile  and  machine  parameters 
shown  in  figure  1.  Note,  the  following  results  are  only  valid  for  those 
times  up  until  the  plasma  first  rebounds  from  its  central  axis.  These 
implosion  times  will  be  given  and  discussed  in  the  results  sections. 


b.  Results  for  variation  of  mass  loading 

In  this  section  we  discuss  the  results  for  varying  the  mass  loading 
while  keeping  the  initial  radial  position  at  1.35  cm.  The  specific  loads 
analyzed  are  5,  10,  20,  30,  40  and  60  /jgms/cm.  Except  for  the  lightest  mass 
loadings,  the  implosion  times  are  long  in  comparison  to  the  time  it  takes 
the  current  to  peak.  Therefore,  for  most  of  the  simulations  discussed  in 
this  section  the  electrical  circuitry  can  be  assumed  to  be  decoupled  from 
the  plasma. 

Figure  2  displays  total,  bound-bound,  bound-free  and  free-free 
radiative  yields.  Similarity,  figure  3  is  almost  identical  to  figure  2, 
except  only  radiation  which  exceeds  1  kev  is  plotted.  From  these  two  plots 
we  find  that  the  components  as  well  as  total  radiation  peaks  at  around  a 
mass  loading  of  20  to  30  }i gms/cm.  In  particular  we  find  that  the  radiative 
yield  of  interest,  i .e.  photon  energy  greater  than  1  kev,  peaks  at  20 
/jgm/cm.  At  this  loading  the  radiative  yield  is  2.2  kilojoules  for  a  plasma 
column  of  4  cms  length.  In  the  remainder  of  this  section  we  want  to 
analyze  figures  2  and  3  in  the  hope  that  we  will  be  able  to  generalize 
these  results  to  other  materials,  mass  loadings  and  circuits. 

AM  of  the  simulations  shown  in  this  section  show  a  (mass)** 
dependence  of  kinetic  energy  upon  mass  loading.  This  result  has  been  noted 
before2,  and  is  due  to  the  fact  that  the  ratio  of  kinetic  energy  stored  in 
a  small  mass  to  that  of  a  larger  mass  is  inversely  proportional  to  their 
masses,  assuming  that  the  force  and  the  time  over  which  the  force  is 
applied  are  equal.  This  stored  up  kinetic  energy  is  transferred  to 
internal  and  radiative  energy  as  the  plasma  assembles  on  axis.  This  is 
readily  seen  in  figure  4  where  peak  values  for  internal,  ionization, 
kinetic,  boundary  work  and  radiative  energies  are  displayed.  The  values 
for  kinetic  energy  are  taken  just  before  the  pinch  whereas  internal  and 
ionization  energies  peak  right  at  the  pinch.  Because  radiative  yield  and 
boundary  work  are  quantities  that  are  integrated  in  time,  their  values  are 
taken  just  after  the  implosion.  Note,  the  drop  in  kinetic  energy  for  the  5 
/jgm/cm  loading  is  because  the  mass  is  so  light  it  implodes  before  the  peak 
in  the  open  circuit  voltage  profile,  see  figure  5. 

Just  because  the  lighter  masses  are  able  to  store  more  kinetic  energy, 
does  not  mean  that  they  will  be  able  to  transfer  this  energy  to  radiation 
at  the  time  of  implosion.  One  reason  this  is  true  is  that  the  (mass)"* 
dependence  also  implies  that  there  will  be  less  mass  to  radiate  during 


22 


implosion.  This  is  especially  true  for  line  radiation  which  is  directly 
dependent  upon  the  number  of  ions  in  an  excited  state.  This  point  is  also 
made  clear  in  figure  6  where  we  see  that  while  peak  power  per  unit  mass 
does  increase  with  smaller  loadings,  the  total  power  eventually  decreases. 
Secondly,  because  the  lighter  masses  are  imploding  with  much  higher 
velocities,  see  figure  7,  they  will  be  spending  much  less  time  near  the 
axis  where  temperatures  and  densities  are  such  that  radiative  power  is 
maximized.  This  is  illustrated  in  figure  5  where  we  see  that  the  "dwell11 
time,  i .e.  the  time  between  peak  values  of  kinetic  energy  before  and  after 
implosion  varies  from  tens  of  nanoseconds  for  the  heavier  masses  to 
fractions  of  a  nanosecond  for  the  lighter  loads. 

We  therefore  find  that  while  it  is  true  that  we  can  couple  more  energy 
into  a  lighter  load  one  eventually  reaches  a  small  enough  loading  that  it 
is  no  longer  possible  to  extract  this  energy  during  assembly  on  axis 
because  of  reduced  mass  and  shorter  assembly  time.  This  point  is  clearly 
illustrated  in  figures  8  and  9  where  we  see  the  almost  total  conversion  of 
kinetic  energy  to  radiation  for  the  30  fi gm/cm  loading  whereas  the  10  pgm/cm 
loading  interaction  is  nearly  elastic.  Note,  how  much  longer  the  dwell 
time  of  the  30  p gm/cm  mass  is  than  the  10  /jgm/cm  mass. 

In  the  beginning  of  this  section  it  is  mentioned  that  the  heavier  mass 
loadings  implode  late  enough  in  time  that  they  can  be  assumed  to  be 
decoupled  from  the  fields  at  the  time  of  implosion.  This  is  seen  in  figure 
8  whereby  if  one  adds  up  the  kinetic  and  internal  energy  just  before 
implosion  and  compares  it  to  the  internal  energy  right  at  implosion,  we 
find  that  the  two  are  approximately  equal.  However,  this  is  not  true  for 
the  10  fi gm/cm  mass  depicted  in  figure  9.  There,  we  discover  that  the 
comb 'nation  of  internal  and  kinetic  energy  before  implosion  is 
significantly  less  than  the  peak  value  of  internal  energy.  The  implication 
is  that  the  lighter  loads  are  imploding  early  enough  in  time  that  there  is 
a  lot  of  magnetic  field  energy  left  that  can  do  boundary  work  on  the  plasma 
during  assembly  on  axis.  Some  of  this  work  will  change  the  internal  energy 
of  the  plasma  and  part  will  be  radiatied  out  of  the  system.  For  the 
example  in  figure  9,  we  find  it  predominantly  changes  the  internal  energy 
of  the  plasma.  The  consequences  of  choosing  one  time  of  implosion  over 
another  will  be  discussed  in  the  next  section. 


Results  for  variation  of  initial  radial  position 
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We  now  fix  the  mass  loading  at  20  figm/cm  and  vary  the  initial  radial 
position  of  the  puff.  The  specific  values  chosen  for  the  following 
simulations  are  centered  at  .25,  .5,  1.0  and  1.35  cms. 

An  examination  of  figures  10  and  11  show  that  we  can  expect  an 
increase  in  the  radiative  yield  of  high  energy  photons,  greater  than  1  kev, 
by  a  factor  of  two  over  the  results  of  the  previous  section  by  changing  the 
initial  radial  position  from  1.35  cms  to  .5  cms.  This  is  an  interesting 
result  because  we  have  the  same  mass  and  dwell  time  but  less  peak  kinetic 
energy  for  the  plasma  starting  at  .5  cms.  The  dwell  time  is  displayed  in 
figure  12  and  the  peak  kinetic  energy  is  shown  in  figure  13.  Therefore,  by 
the  arguments  of  the  previous  section,  we  expect  figure  13  to  show  that 
radiative  yield  peaks  at  the  same  radial  position  as  kinetic  energy,  i .e. 
at  an  initial  position  of  1  cm.  However,  what  we  find  instead  is  that 
radiative  yield  is  maximized  by  starting  even  closer  to  the  axis  at  an 
initial  radial  position  of  .5  cm. 

The  discrepancy  lies  in  the  fact  that  the  arguments  of  the  last 
section  are  based  on  the  premise  that  the  plasma  is  decoupled  or  nearly 
decoupled  from  the  electrical  fields  during  assembly  on  axis.  We  can  see 
from  the  implosion  times,  see  figure  12,  that  the  only  plasmas  nearly 
decoupled  from  the  field?  r»<'e  those  with  inital  positions  greater  than  1 
cm.  There  is  a  different  phenomenon  taking  place  for  plasmas  with  smaller 
initial  positions  because  they  are  still  strongly  coupled  to  the  fields. 
What  is  happening  is  that  the  decoupled  plasmas  encounter  only  a  single 
shock  as  they  hit  the  axis  whereas  the  coupled  plasmas  are  multiply 
shocked.  These  coupled  plasmas  are  not  only  shocked  by  running  into  the 
axis  but  they  continue  to  be  shocked  by  the  fields  as  boundary  work  is  done 
on  the  plasma.  The  situation  is  similar  to  inertial  confinement  fusion 
(ICF)3  where  it  is  well  known  that  if  you  want  to  compress  a  fluid  it 
should  be  done  by  multiple  shocking  along  an  isentrope  rather  than  one  big 
shock.  The  lower  the  adiabat  of  the  shock  the  larger  the  compression  for 
the  same  imploding  pressure.  This  is  the  reason  that  the  coupled  plasmas 
having  lower  peak  kinetic  energy  have  correspondingly  lower  temperatures 
and  higher  densities  at  implosion,  see  figure  14.  The  lower  kinetic  energy 
means  less  preheating,  therefore  the  shock  adiabat  is  lower. 

The  way  this  relates  to  an  increase  in  radiation  for  the  smaller 
radii  plasmas  is  that  recombination  and  bremsstrah I ung  radiation  are 
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proportional  to  the  square  of  the  density.  We  therefore  expect  the 
radiative  yield  to  increase  for  plasmas  that  are  more  strongly  coupled  to 
the  fields  because  they  are  more  compressed.  The  plasma  will  be  most 
strongly  coupled  if  it  implodes  right  about  the  time  of  maximum  voltage 
which  for  these  simulations  corresponds  to  the  .5  cm  case.  Although  the  .25 
cm  simulation  is  imploding  so  early  in  time  that  the  fields  have  not  had 
time  to  gain  maximum  strength.  The  kinetic  energy  is  low  enough,  see  figure 
13,  that  the  preheating  it  produces  is  small  compared  to  the  other 
simulations.  This  allows  for  compression  to  much  higher  densities  at  the 
expense  of  lower  temperatures  as  is  illustrated  in  figure  14.  The  point  is 
also  made  in  figure  10  where  we  see  that  there  is  only  a  difference  of  600 
joules/cm  of  total  radiative  yield  between  the  .25  cm  and  .5  cm 
simulations,  but  an  examination  of  figure  11  shows  that  all  of  it  comes 
from  the  high  energy  component,  corresponding  to  higher  temperatures. 
Another  illustration  of  the  density  dependence  of  radiation  is  again 
depicted  in  figures  10  and  11  where  one  sees  a  dramatic  increase  in 
recombination  radiation  but  a  drop  in  line  radiation  as  the  initial  plasma 
radius  is  moved  inward  allowing  the  plasma  to  be  more  strongly  coupled  to 
the  fields. 

d.  Summary  and  conclusions 

We  have  shown  that  imploding  gas  puffs  operate  in  two  regimes,  one  in 
which  the  plasma  implodes  late  enough  in  time  that  it  can  be  considered  to 
be  decoupled  from  the  electrical  fields  as  it  hits  the  axis  and  another  in 
which  the  plasma  implodes  early  in  time  such  that  it  is  still  coupled  to 
the  fields. 

For  the  plasmas  that  are  decoupled  we  find  that  more  kinetic  energy 
can  be  stored  in  a  lighter  mass.  However  this  energy  is  less  efficiently 
converted  to  radiation  not  only  because  of  the  smaller  mass  but  also 
because  the  plasma  spends  less  time  in  its  maximum  imploded  state.  If  one 
is  trying  to  maximize  the  radiative  yield  for  a  plasma  that  is  decoupled 
from  the  electrical  fields  it  is  a  trade  off  of  getting  enough  energy  into 
the  system  and  then  being  able  to  efficiently  radiate  at  implosion.  For  an 
argon  plasma  that  starts  sufficiently  far  away  from  the  axis  that  it  is 
decoupled  from  the  fields  at  implosion,  we  found  that  the  best  radiative 
yields  where  achieved  using  a  mass  loading  of  20  to  30  ji gms/cm.  At  this 


setting  we  obtain  total  radiative  yields  of  4.8  kilojoules  of  which  almost 
half  is  composed  of  photons  with  energies  greater  than  1  kev . 

The  coupled  plasmas  exhibit  a  different  behavior  that  is  similar  to 
ICF  in  that  the  fields  multiply  shock  compress  the  plasma  as  it  implodes. 
The  adiabat  upon  which  the  plasma  is  compressed  is  dependent  upon  the 
kinetic  energy  of  the  plasma,  the  lower  the  kinetic  energy  the  lower  the 
adiabat  and  consequently  the  compression  is  higher.  Because  the  plasmas 
starting  closer  to  the  axis  have  less  time  to  gain  kinetic  energy,  they 
become  more  compressed.  Since  recombination  and  bremsstrahlung  radiation 
are  proportional  to  the  sqare  of  the  density,  the  radiative  yield  increases 
for  plasmas  starting  closer  to  the  axis  because  of  their  higher  densities. 
The  adiabat  upon  which  the  plasma  is  compressed  determines  the  radiative 
spectrum.  We  find  for  argon  that  the  radiative  yields,  both  total  and  high 
energy,  are  doubled  over  the  results  previously  mentioned  if  the  initial 
starting  position  of  the  plasma  is  such  that  the  plasma  implodes  on  axis 
about  the  same  time  as  the  open  circuit  voltage  profile  peaks. 

Note,  more  investigation  needs  to  be  done  on  the  plasmas  that  are  coupled 
to  the  fields  because  they  are  in  a  regime  that  is  susceptible  to 
instab i I ities. 

All  of  the  simulations  presented  in  this  work  are  for  a  Gamble  II  - 
like  open  circuit  voltage  profile  that  corresponds  to  the  absence  of  a 
plasma  erosion  opening  switch  (PEOS) .  Since  it  has  been  experimentally 
found  that  the  presence  of  a  switch  lowers  the  peak  current,  the  arguments 
of  this  report  predict  that  the  peak  velocities  and  radiative  yields,  for 
a  given  peak  current,  would  also  be  reduced.  This  prediction  is  in 
agreement  with  "SIMPLODE"  ^  a  0-D  radiation  hydrodynamic  model  developed 
by  the  Plasma  Radiation  Branch  that  simulates  both  the  switch  and  non¬ 
switch  cases.  However  this  prediction  is  only  valid  when  the  fields  are 
decoupled  from  the  plasma  at  implosion,  more  analysis  needs  to  be  done  for 
the  coupled  plasmas.  Also,  this  prediction  does  not  include  any  stability 
analysis,  which  needs  to  be  included  in  order  to  make  a  judgement  about  the 
merits  of  a  switch  versus  a  non-switch.  It  is  conceivable  that  the 
benefits  of  increased  uniformity  of  the  plasma  resulting  from  faster  rise 
times  produced  by  the  use  of  a  PEOS  far  outway  the  side  effects  of  reducing 
the  plasma  current. 


Comparisons  of  the  results  of  this  model  with  SIMPLOOE  are  in 
excel lant  agreement  when  it  comes  to  global  comparisons  of  radiative  yield 
versus  mass  loading  for  a  fixed  initial  radius,  i .e.  assuming  that  the 
fields  are  decoupled  at  implosion.  Although  the  SIMPLODE  results  are  for 
neon  the  trends  are  almost  identical.  It  is  therefore  recommended  that 
that  the  0-D  models  be  utilized  for  overall  scoping  of  a  problem  and  then 
for  finer  details,  such  as  producing  a  radiative  spectrum,  the  1-D  models 
should  be  employed. 
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Fig.  1  —  Gamble  II  voltage  profile  and  circuit  diagram. 
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Fig.  7  —  Peak  A  -  electron  temperature,  B  -  ion  density,  and  C  - 
velocity  versus  mass  loading,  temperatures  and  densities  peak  at 
implosion  while  velocities  peak  before  implosion. 
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Fig.  10  —  Radiative  yield  versus  initial  radial  position,  A  is  total,  B  is 
bound-bound,  C  is  recombination  and  D  is  bremsstrahlung  radiation. 
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Fig.  1 1  —  K-shell  radiative  yield  versus  initial  radial  position.  A  is 
total.  B  is  bound-bound,  C  is  recombination  and  D  is  bremsstrahlung 
radiation. 
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Fig.  13  —  Peak  values  for  A  -  internal.  B  -  ionization.  C  -  kinetic.  D  - 
boundary  work  and  E  -  radiative  yield  energies  versus  initial  radial 
position. 


Importance  to  Uniformity  in  a  Z-Pinch 


In  order  for  a  Z-pinch  to  exhibit  a  narrow  radiation  signature  in  time, 
the  plasma  column  must  be  uniform  along  the  Z  direction  during  the  time 
when  the  optimum  implosion  conditions  are  realized.  Of  particluar 
importance,  is  the  manner  in  which  the  hollow  cylindrical  shell  of  gas 
formed  by  a  high  Mach  number  nozzle  behaves  during  the  runin  phase  of 
the  implosion.  The  spread  in  the  gas  due  to  expansion  after  exiting  the 
nozzle  can  make  the  stagnation  time  of  the  implosion  nonuniform  in  Z. 
This  phenomena  is  often  referred  to  as  "zippering”.  In  addition  to  the 
annular  outside  puff,  these  configurations  may  also  include  a  central  gas 
puff  which  can  serve  as  a  radiation  convertor1.  When  the  outside  puff 
stagnates  on  the  inner  puff,  the  zippering  effect  is  transferred  to  the  inner 
puff  resulting  in  density  variations  along  Z  at  the  time  of  peak 
compression.  In  this  paper,  we  will  examine  representative  nozzle  designs 
in  use  on  experiments  which  are  currently  being  planned  and/or  conducted 
at  the  Physics  International  Company  (  hereafter,  referred  to  as  PI  ). 
Both  the  outside  and  central  puff  are  taken  to  be  krypton.  The  mass  of 
the  inside  puff  is  set  to  a  variable  fraction  of  the  outside  puff  mass.  The 
three  designs  which  have  been  considered  are: 

(1)  Outside  radius  of  2.2  cm,  inside  radius  of  0.8  cm,  and  outside  puff 
centered  at  a  radius  of  1.4  cm.  Angular  spread  of  20  degrees  on  the 
outside  and  angular  spread  of  9  degrees  on  the  inside  (  i.e.,  inside  meaning 
on  the  side  closest  to  the  center  of  the  column  ).  The  differences  in 


angular  spread  have  been  attributed  by  PI  to  interactions  with  the  central 
gas  puff.  This  gives  rise  to  less  of  a  spread  on  the  inside. 

(2)  Outside  radius  of  2.2  cm,  inside  radius  of  1.8  cm,  and  centered  at  2.0 
cm.  For  the  outer  puff,  an  outside  spread  of  25  degrees  and  an  inside 
spread  of  15  degrees.  Expansion  of  the  inside  puff  is  20  deg  ees. 

(3)  Same  as  (2)  above  but  with  an  overall  tilt  of  the  outside  puff  inward. 

The  radius  of  the  inner  puff  is  always  taken  to  be  2.5  mm.  Unless 
otherwise  stated,  the  inner  puff  can  be  assumed  to  have  an  angular 
expansion  of  20  degrees.  Figure  1  shows  a  schematic  of  the  nozzle  design 
and  further  illustrates  the  meaning  of  these  values.  In  this  report,  we  will 
specifically  discuss  designs  (2)  and  (3). 

Computational  Considerations 

These  simulations  were  carried  out  using  a  t\  o-dimensional  Lagrangian 

hydrodynamics  code.  Earlier  applications  and  discussions  of  this  code, 

2 

PRISM,  appear  in  previous  annual  reports  .  The  one-dimensional  version 
of  this  code  was  used  to  investigate  the  radiative  yields  of  double  puff 
implosions  and  was  reported  on  in  reference  1. 

The  typical  mesh  size  was  30x40  with .  spacing  determined  by  resolution 
considerations  in  the  two  puffs.  The  nozzle  flow  itself  is  neglected  except 
in  the  initial  density  profiles.  The  density  variations  are  given  by  Gaussian 
profiles  in  the  radial  direction  with  the  full-width-at-half-maximum 
increasing  along  Z  according  to  the  angular  spreads  dictated  by  the 


particular  nozzle  design.  The  angles  chosen  are  described  above.  After 
computing  the  profiles  in  both  R  and  Z,  the  density  is  normalized  to  a 
given  mass  per  unit  length  for  each  slice  along  Z.  The  mesh  is  extended 
at  both  the  top  and  bottom  of  the  pinch  in  order  to  prevent  disturbances 
which  propagate  in  Z  from  being  reflected  back  into  the  pinch  region  (  i.e., 
Z  between  0  cm  and  3  cm  ). 

In  addition  to  the  hydrodynamics,  two-temperature  thermal  transport 
with  an  idealized  equation  of  state  is  used.  The  magnetic  field  is  advanced 
according  to  resistive  diffusion  and  a  boundary  value  given  by  a  temporal 
current  profile. 


Discussion  of  the  Results 


Several  computer  runs  were  made  in  which  the  initial  nozzle  parameters 
as  well  as  the  ratio  of  inside  to  outside  mass  were  varied.  The  parameters 
are  those  of  designs  (2)  and  (3),  which  where  discussed  above.  Both  the 
temporal  profile  and  magnitude  of  the  current  were  fixed.  They  were 
provided  for  these  simulations  by  PI.  Our  main  interest  in  these  runs  was 
to  look  for  regions  of  uniform  density  along  or  near  the  central  axis.  For 
this  reason,  a  graphics  contour  plotting  routine  which  displayed  only  the 
50%  and  25%  values  of  the  peak  electron  density  was  used.  The  results 
from  the  hydrodynamic  simulations  are  discussed  below. 

Figure  2(a)  shows  the  density  contours  near  the  peak  of  the  implosion 
of  a  double-puff  design  in  which  the  central  puff  had  a  mass  equal  to  80% 
of  the  outside  mass.  In  this  run,  the  tilt  angle  was  set  at  5  degrees.  The 
50%  of  peak  density  contour  is  shown  with  a  solid  line,  while  the  25% 
contour  is  dashed.  The  inside  contour  is  seen  to  have  about  1  cm  along 
the  central  axis.  The  length  of  the  25%  contour,  however,  is  1.5  cm. 


Both  contours  are  contained  within  a  1  mm  radius  (  note  that  the  R  ana 
Z  axis  scales  are  not  the  same  ). 

Figure  2(b)  shows  a  case  in  which  the  inner  puff  contains  20%  of  the 
mass  as  that  of  the  outer  puff.  In  this  example,  there  is  no  tilt  to  the 
outer  puff.  It  can  be  seen  that  the  length  of  the  50%  contour  is  roughly 
1.5  cm  while  the  25%  contour  is  over  2  cm  in  length.  The  central  portion 
of  the  inner  contour  is  also  centered  lower  in  Z  than  the  case  discussed 
above.  It  has  been  seen  in  these  runs  that  when  the  expansion  angle  of 
20  degrees  is  used  in  the  inside  puff,  the  tilted  outside  puffs  tend  to  be 
better  centered  in  Z.  However,  the  overall  length  along  Z  does  not 
necessarily  correlate  with  tilt. 

In  contrast,  Fig.  2(c)  shows  a  case  in  which  there  is  no  expansion  of 
the  inner  puff.  Expansion  of  the  inner  puff  may  be  controllable  with 
proper  nozzle  design.  The  use  of  no  expansion  represents  a  somewhat 
idealized  case.  All  other  parameters  are  the  same  as  those  of  Fig.  2(a). 
Without  the  expansion,  the  central  puff  is  more  elongated  and  uniform 
along  Z  providing  for  a  more  uniform  density  distribution.  In  fact,  the 
length  of  the  inner  contour  extends  well  over  2  cm  along  the  center  of  the 
pinch. 

Finally,  Fig.  2(d)  shows  temperature  contours  corresponding  to  the 
example  shown  in  Fig.  2(b).  In  this  particular  instance,  temperatures  of 
only  300-500  eV  were  reached  during  the  implosion  in  the  regions  shown  in 
Fig.  2(b).  In  other  cases,  temperatures  as  high  as  800  eV  have  been  seen 
but  not  under  the  optimal  density  conditions  displayed  in  Fig.  2.  In  any 
case,  the  temperature  may  be  somewhat  suspect  because  of  the  crude 
ionization  and  thermal  energy  partition  employed  in  these  2D  runs.  What 
can  be  seen  in  the  plots  is  that  the  temperature  gradient  is  at  least  as 
smooth  as  the  density  gradient.  This  implies  that  the  radiative  output 


& 


should  be  nearly  uniform  in  this  region.  Whether  this  produces  increased 
yields  of  hard  K-shell  radiation  is  yet  to  be  determined.  It  may  increase 
the  total  yield,  however. 


Conclusion 


The  hydrodynamic  simulations  presented  here  indicate  that  regions  of 
relative  uniformity  can  be  produced  in  Z-pinch  implosions.  Lengths  of  up 
to  2  cm  in  the  Z  direction  at  a  density  equal  to  50%  of  the  peak 
implosion  density  have  been  demonstrated  within  the  context  of  these 
simulations.  Temperature  gradients  are  also  seen  to  have  long  scale  lengths 
in  these  regions. 

Although  radiative  yields  were  not  themselves  investigated  in  this 
study,  the  implication  is  that  certain  nozzle  designs  can  produce  large 
spatial  regions  at  implosion  which  would  lead  to  radiative  outputs  more 
peaked  in  time.  Hence,  higher  radiative  power  levels  can  be  produced  for 
a  given  total  yield.  It  is  also  probable  that  greater  total  yields  will  result 
from  the  more  uniform  implosions. 
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Figure  1.  Schematic  of  the  double-puff  nozzle  design.  (a)  radius  of  the 

inner  puff,  (b)  inside  radius  of  the  outer  puff,  and  (c)  outside 
radius  of  the  outer  puff.  The  angular  expansions  are  shown  by 
the  lines  drawn  from  points  a,  b,  and  c  to  the  anode  (right 
hand  side  of  the  figure).  In  this  figure  no  tilt  is  shown  for 
the  outside  puff. 
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Section  E.  Opening  Switch  and  Diode  Modeling 


Gyrokinetic  Propagation  Studies 

The  behavior  of  collsionless  peripheral  plasma  around  PRS  loads  or  in  the  highly 
magnetized  regions  of  VOS  devices  represents  a  crucial  area  in  one’s  physical  understanding 
of  the  power  flow  in  such  devices.  In  the  radial  transmission  line,  particularly,  the  existence 
of  a  tensor  pressure,  and  the  rapid  radial  acceleration  and  compression  of  guiding  centers 
forces  a  very  complex  ”Ohms  Law”  on  the  plasma  medium  in  the  collisionless  limit.  As  a 
prelude  to  examining  the  detailed  microstability  of  such  plasmas,  a  one  dimensional  radial 
propagation  study  has  been  undertaken  in  1986. 

The  propagation  calculation  has  two  major  goals: 

(i)  to  assess  the  transmission,  reflection,  and  absorption  properties  of  collisionless 
gyrokinetic  plasmas  in  a  radial  transmission  line,  and, 

(ii)  to  determine  the  situations  where  microinstabilities  are  likely  modes  of  plasma 
turbulence  in  the  tenuous  regions. 

The  propagation  computation  is  summarized  below  and  in  prior  work.  Owing  to  the 
virulently  nonlinear  character  of  the  current  density  in  the  gyrokinetic  limit,  the  studies 
required  to  reach  the  first  of  these  goals  are  still  ongoing.  However  those  results  fall  into 
place,  the  qualitative  behavior  of  the  field  profiles  which  thereby  emerge  has  begun  to  shed 
some  light  on  the  second  objective. 

This  calculation  rests  on  four  major  elements: 

(i)  a  smooth  set  of  radial  mass  density,  magnetic  and  electric  field  profiles;  consistent 

with 

(ii)  a  profile  of  mean  perpendicular  kinetic  energy  density  determined  by  a  Fokker- 
F’lanck  calculation  of  the  ’’driftaway”  process; 

(iii)  an  exterior  and  interior  circuit  model  delivering  energy  to  the  annular  gyrokinetic 
plasma;  and 

( i v )  a  fully  implicit,  fully  electromagnetic.  1-1)  Lagrangian  method  of  lines  solver  for 
the  field  and  fluid  equations  which  advances  the  guiding  center  drift  speed  and  the  fluid 
Jacobian  equivalent  to  Ampere's  and  Faraday’s  Laws  in  the  collisionless  gyrokinetic 
system. 

Kach  of  these  elements  will  he  discussed  and  some  general,  but  preliminary  conclusions 
reported. 

Exact  separable  implosion  with  n  /  />  . 


For  special  initial  conditions  there  exist  exact  solutions  to  the  gyrokinetic  moment 
equations  that  describe  strongly  magnetized  plasma  in  the  periphery  of  a  z-pinch.  An 
isotropic  pressure  implies  that  the  collisions  are  sufficient  to  keep  the  velocity  components 
randomized.  This,  then,  implies  a  finite  conductivity,  in  contrast  to  the  assumption  of 
infinite  conductivity  which  demands  a  collisionless  plasma.  However,  as  the  plasma  gets 
closer  to  the  axis  in  a  cylindrical  geometry  the  pressure  pj|  parallel  to  the  magnetic  field 
B  =  Bg(r)  increases  as  l/r2,  while  the  perpendicular  pressure  increases  as  1  /r.  In  fact, 
in  any  curvilinear  magnetic  field  an  isotropic  pressure  is  strictly  speaking  inconsistent  with 
infinite  conductivity. 

The  fluid  equation  of  motion  in  the  radial  direction,  for  an  anisotropic  pressure,  in¬ 
cludes  an  centrifugal  force  Z,2/r3  that  reflects  the  conservation  of  angular  momentum, 


DU  pL  n  1  d{rpBWL) 

P  n.  —  3  ~~  ^  j zB 

Dt  rJ  r  dr 

The  macroscopic  quantities  that  enter  into  the  moment  equations  are  averages  over 
the  particle  distribution  function.  Without  collisions  the  distributions  are  not  necessarily 
maxwellian,  but  the  exact  functional  form  is  unknown.  Instead  of  the  distribution  in 
terms  of  the  velocity  components  in  cylindrical  geometry  /(r, v,t)  it  is  better  to  use  the 
distribution  function  in  terms  of  the  dynamical  invariants,  the  magnetic  moment  p  and 
the  angular  momentum  P$.  The  two  equivalent  forms  of  the  distribution  function  are 
connected  by 


/(r,  v,  <)  = 


m2  r 


n{r,t)  F{p,Pg,r)  . 


[2tt  B(r,t)  J 

The  factor  in  '•  •  ■!  is  the  Jacobian  in  the  velocity  space  transformation  to  p  and  Pg.  Without 
collisions  the  invariants  for  each  particle  do  not  change,  and  the  distribution  function  F  is 
constant  in  time.  The  particle  density  n(r,  t)  is  no  longer  part  of  F,  which  is  normalised  as 


m2  r 


27r  B(r,t) 


J  dpdPg  F(p,  Pg,r)  =  1, 


for  both  electrons  and  ions.  With  the  assumption  of  quasineutrality  the  mass  density  p  is 
p  -  ;  mt„n,,)n(r, /). 

The  parallel  pressure  is 


j  dpd  Pg  F P#  j  m2, 


and  the  perpendicular  pressure  becomes 


pH', 


J  dpd  Pfi  Fp/rn , 


with  a  summation  over  species  implied. 


Equation  (1)  with  vanishing  angular  momentum,  L\\  =  0,  yields1  a  selfsimilar  solution, 
but  it  turns  out  that  even  for  non-zero  L||  a  whole  class  of  seif-similar  solutions  are  possible. 
The  one  highlighted  here  is  relevant  for  electromagnetic  propagation  modeling. 

Henceforth  the  dependent  variables  B,p,  L\\,  W±  and  their  derivatives  are  functions  of 
the  single  selfsimilar  variable  £  =  r/a(t),  where  a(t )  is  a  typical  scale  length  that  contains 
the  sole  time  dependence.  Then  Eq.  (1)  becomes 

d2a  =  Hi  _  wj 

dt2  a 3  a2  a  ’ 

where  a(t)  is  the  normalized  scale  length  a(t)  =  a(t)/ao. 

Here  w||2(£)  =  T||(0/£4  depends  in  general  on  the  spatial  variable  £,  and  Eq.  (2) 
does  not  separate.  Separation  is  feasible  only  when  the  initial  angular  momentum  is 
proportional  to  (radius)4,  i.e.  w2  is  constant.  Recall  that  L 2  is  an  average  of  P2,  where 
Po  ~  v$r  is  the  particle  angular  momentum.  Since  v$  ~  u$r  also  P$  ~  ui$r2:  constant  ui2 
then  implies  a  parallel  temperature  oc  r2. 

The  perpendicular  pressure  term  gives  rise  to 

2  _  1  d\Hs.mp] 

which  should  no  longer  be  a  function  of  £.  In  this  term  the  spatial  dependence  of  the  density 
p{£)  can  be  played  off  against  the  perpendicular  thermal  energy  H±{£)  =  5(£)W_l(£)- 
Either  p  or  H_ i  can  be  chosen  freely.  It  is  easiest  to  separate  the  spatial  dependence  in 
//jl  by  setting  H±(£)  =  Ho/g{£)-  The  mass  density  is  then 


The  acceleration  term  contains  the  magnetic  field  through  u>a,  which  is  the  transit 
frequency  of  an  Alfven  wave  though  the  pinch’  scale  size  a 0,  u>\  =  B2 /4npa2.  A  constant 
then  implies  a  specific  magnetic  field,  in  terms  of  g(£): 

=  exp  dxx9{x)  , 

where  B ^  =  47ra2w^Poi  the  magnetic  field  at  the  characteristic  distance  a0.  This  equation 
also  defines  the  total  current  enclosed  by  the  pinch  within  a  distance  f,  because  ££  = 
2 I($)/ca0. 


The  function  <7(£)  is  arbitrary  within  restrictions  defined  by  certain  physical  require¬ 
ments.  One  is  that  the  density  p  be  positive  and  finite  everywhere,  another  is  that  the 
various  w2’s  be  positive.  The  choice  here  is  g(£)  =  £2,  which  satisfies  all  these  require¬ 
ments.  With  this  special  form  for  g(£)  the  various  integrations  give  rise  to  known  functions. 
For  example,  the  density  becomes 

p  =  p0el/4£exp-  £4/4, 

where  p0  is  given  by  the  mass  per  unit  length  within  £,  A(£),  through 

A  =  2tt  alp0Wl/A{£), 

where  is  the  incomplete  gamma-function  of  order  1/4.  Likewise,  the  current 

carried  by  the  corona  enclosed  within  £  becomes 

I2{t)  -  /2(0)  =  fl^e1/4  f  dxx4exp(-x4/ 4), 

Jo 

where  7(0)  represents  the  current  carried  by  a  collisional  pinch  located  on  axis. 

In  any  particular  circuit  characterized  by  the  current  risetime  and  the  interior  resis¬ 
tance  of  the  bulk  plasma  load,  one  can  use  the  results  of  the  Fokker-Planck  calculation 
to  establish  a  scale  sound  speed  H0  for  the  selfsimilar  startup  profiles,  using  the  initial 
voltage  and  normalized  magnetic  field  fi  implied  in  any  circuit  model. 

Shown  in  Figure  1  is  the  variation  of  H0  =  H x/ H with  the  parame¬ 
ter  [E  /  E  p)  /  [Ut  0)  at  30  t0  (collision  times)  for  some  190  combinations  of  E/Ep  and 
Qt0  .  The  clearly  one  dimensional  variation  of  the  Fokker-Planck  result  is  not  manifest 
a  priori,  but  it  certainly  helps  reduce  the  effort  in  scanning  initial  conditions.  Shown  in 
Figure  2  is  an  initial  condition  profile  derived  from  these  considerations. 

Collisionless  current  density. 

In  a  collisional  plasma  the  collisions  are  the  dominant  influence  that  balances  the 
motion  of  the  plasma  electrons  against  the  acceleration  by  the  electric  field.  The  resulting 
collisional  electron  drift,  and  thus  the  current  density  j,  is  proportional  to  the  electric 
field  E,  j  =  aE.  The  electrical  power  input  into  the  plasma  ends  up  in  increased  random 
velocities  of  the  plasma  particles,  Ohmic  heating. 

Even  in  the  absence  of  collisions  the  magnetized  plasma  in  the  z-pinch’  periphery  can 
absorb  energy  from  the  electrical  pulse.  In  this  case  the  energy  does  not  go  into  random 
motion  of  the  plasma  particles,  but  instead  is  put  into  organized  motion.  A  strongly 
magnetized  plasma  in  an  electric  field  is  not  accelerated  by  the  field,  but  drifts  with  a 
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velocity  v  =  E / B.  For  constant  E  the  energy  ’-emains  constant,  but  when  the  electric  held 
increases  in  time  the  plasma  accelerates.  The  kinetic  energy  density  pv2  increases  also  (p 
is  the  plasma’s  mass  density),  and  electrical  power  is  absorbed.  Likewise,  electrons  that 
E  \  B  drift  into  a  region  of  increased  magnetic  field  spin  up,  increasing  their  perpendicular 
energy  in  their  effort  to  keep  the  magnetic  moment  p  ~  w L/ B  constant.  The  increase 
in  w i,  or  the  increased  magnetization,  is  another  power  sink.  In  a  cylindrically  symmetric 
geometry  with  B  ~  1  / r  the  particle's  angular  momentum  about  the  axis,  P$  =  mrvs,  is 
conserved:  a  drift  toward  the  axis  increases  the  kinetic  energy  parallel  to  the  magnetic 
field  u>||  as  1  /r2 . 

The  axial  current  density,  calculated  either  by  momentum  balance  or  by  averaging  of 
net’*,  is  then 

pc2  D{E  B)  c 
1  B  I)t  '  B 

The  first  term  is  related  to  the  acceleration  of  the  /v  >  B  drift  motion,  the  second  comes 
from  angular  momentum  conservation,  or  />  .  and  the  third  is  the  magnetization  con¬ 
tribution,  from  p  In  the  acceleration  term  the  ions  dominate;  the  electrons  are  more 
important  for  the  two  pressure-derived  terms,  particularly  in  the  general  case  without 
thermal  equilibrium  between  electrons  and  ions. 


PL*  \d[rpB\\\) 
r!  r  . dr 


Electromagnetic  energy  transport  through  the  corona. 

What  fraction  of  the  input  power  pulse  passes  through  the  collisionless  but  energy- 
ab:  orbing  plasma  corona?  The  transmission  coefficient  'I'  is  the  power  reaching  the  central, 
collisional  load  at  some  inner  boundary.  S,  2trr,  (E;  *  Bj),  and  the  power  coming  in 
through  the  outer  boundary,  Sh,  i.e. 

e  Et  Hi  ri 

1  Sh  ~  E-  »  : 

Eb  Bb  rb 

the  absorption  coefficient  is  .4  -  1  T.  While  the  inner  boundary  is  nearly  stationary, 
the  outer  boundary  is  moving  with  the  plasma,  r ^  -  rt,(f).  Motion  of  the  plasma’s  outer 
boundary  changes  the  volume  of  the  vacuum  between  the  plasma  and  the  power  inlet.  The 
magnetic  field  energy  in  this  vacuum  changes  correspondingly.  The  transmission  efficiency 
of  the  coroi  a  plasma  including  the  outside  vacuum  is  then  defined  as  t  =  5,/S,  where  S 
is  the  total  power  transmitted  through  the  inlet. 

Toward  the  inlet  the  corona  plasma  goes  over  smoothly  into  hard  vacuum,  while  on 
the  other  side  the  plasma  becomes  denser  until  the  collisions  begin  to  dominate.  This  is 
where  the  central  pinch  starts.  The  inductive  approximation.  dEjdt  47rj,  is  invalid  in 
the  outer  region  without  plasma;  a  proper  transition  between  regimes  then  demands  that 
this  term  be  kept  also  in  the  corona  regime.  The  current  density  in  Eq.  (3)  implies  that 
the  magnetic  field  moves  with  the  plasma,  suggesting  that  a  Lagrangian  coordinate  system 
fixed  in  the  fluid  may  simplify  the  computations. 


The  Lagrangian  radial  coordinate  R  =  R(£,t)  evolves  from  the  coordinate  £  at  time 
t  —  0  according  to  the  implicit  relation  dR/dr  =  —U(R,t),  and  R  is 


R{t,r)  =  Z-[  dT'U{R,T'), 

Jo 


A  fluid  element  d£  changes  according  to  the  Jacobian  J  =  dR/d£,  where  DJ/Dt  = 
-dU{R,r)/dt. 

Magnetic  flux  conservation  requires  /  B£d£, 

B{R,  t)  =  B{£)d£/dR, 


«,V 

•V 

♦»&< 


IWlXJ 


The  mass  in  a  shell  ( pRdR )  is  conserved,  whence 

Rp{R,r)  =  tp(i)d£ldR. 

From  these  relations  the  Alfven  speed  —  B2 /4np  transforms  as 

c\(R,r)=c2A(Od(ln()/d(lnR). 

In  terms  of  the  normalized  plasma  velocity  U  =  E/B  Ampere’s  law  becomes 
DU  B  as(l  -  IP)  1  2 

IT,  -  [r  +  M - 

Here  K\  -  1  +  c2 / c2A(R ,  r)  determines  the  wave  speed  c/if,*:  without  plasma,  c,*  — >  oo, 
and  A' .4  becomes  unity,  for  a  wave  speed  c.  KA  becomes  c/cA  in  the  opposite  case,  the 
wave  speed  is  eA,  and  the  inductive  approximation  is  valid. 

The  current  density  Jr  is  the  part  of  the  current  density  in  Eq.  (3)  coming  from  the 
particle  pressure,  i.e.  the  current  density  without  the  contribution  from  accelerating  the 
plasma  drift: 

B  VF',|  BWi  d[pBW±(d(/dR)}2  £  dR 
Jr  R3  c2a  H'h  ipBWL  R  d£ ' 

The  particle  pressure-related  terms  VF||(£)  and  are  determined  by  the  selfsimilar 

solutions  at  r  =  0  (cf.Fig  2),  and  these  are  used  as  initial  conditions. 

The  exterior  boundary  condition  comes  from  the  electromagnetic  wave  driving  the 
plasma.  Typically  the  electric  field  is  connected  to  some  pulse  power  device,  cf.  Figure  3, 
which  should  in  principle  be  influenced  by  the  load.  The  interior  boundary  condition  is 
determined  completely  by  the  axial  resistance,  viz. 

V  =  R  RLoadjl 

<  <  2Znt„ 


while  the  exterior  condition  is  posed  as  an  alternative  acceleration  formula, 

K\  p  t  =  f{Vs,Lr>iode,Ls,R>,ZD), 

with  the  circuit  current  Is  and  voltage  V"s  determined  in  the  usual  way.  The  derivation 
of  this  circuit  equation  presumes  continuity  in  the  line  impedance  at  the  plasma  vacuum 
interface;  the  next  level  of  approximation  is  to  relax  this  constraint  and  go  to  a  transmission 
line  model.  The  capacitance  and  inductance  of  the  circuit  model  then  set  the  circuit 
risetime;  a  reasonable  risetime  domain  is  t,,  —  (10,50)  ns. 

Shown  in  Figure  4  is  a  typical  time  development  of  the  plasma  electric  field,  magnetic 
field,  and  current  density  as  a  50  ns  risetime  pulse  is  impressed  on  the  gyrokinetic 
corona.  The  energy  is  not  absorbed  gracefully;  the  plasma  current  density  in  Fig.  4c  shows 
wiggles,  which  are  TEM  waves  (apparently  analogous  to  magnetosonic  MHD  modes)  that 
are  excited  by  the  relatively  sudden  application  of  the  electric  field.  For  low  values  of  H0 
the  profiles  change  little  over  time,  and  the  electric  field  is  easily  transmitted.  For  larger 
values. e.g.  Hp  =  (7.0  •  10Cl)2  — ►  Te  -  20eV  ,  the  electric  field  is  heavily  modified  and  the 
input  power  is  substantially  attenuated  by  the  corona. 

Summary  and  conclusions. 

Current  can  be  lost  in  the  periphery  of  a  z-pinch  by  various  mechanisms.  Initially  the 
z-pinch  carries  no  current,  and  the  inductive  voltage  at  the  edge  of  the  cathode  away  from 
the  pinch  material  can  drive  a  currrent  in  free  electrons.  The  free  electrons  can  be  the 
bulk  of  the  current  at  the  very  start,  but  this  current  is  cut  off  rapidly  once  the  current 
starts  flowing  through  the  z-pinch. 

The  tenuous  plasma  outside  the  z-pinch  can  also  carry  current.  As  a  remnant  of  back¬ 
ground  gas  or  perhaps  the  edge  of  the  injected  gas  column  in  a  puff-gas  z-pinch,  this  plasma 
is  not  well-specified  either  experimentally  or  theoretically.  How  much  current  is  shunted 
through  this  plasma  depends  linearly  on  the  plasma  density,  and  also  on  parameters  such 
as  the  temperature  and  the  magnetic  field  in  the  periphery,  i.e.  the  current  in  the  main 
z-pinch. 

The  plasma  temperature  in  the  periphery  at  the  end  of  "drift away”  is  computed  from 
a  simplified  Fokker-Planck  analysis,  assuming  that  the  electron  velocity  remains  isotropic 
through  collisions  with  the  ions.  Plasma  heating  proceeds  until  the  electron  collsion  time 
is  comparable  to  the  pulse  time,  resulting  in  an  upper  temperature  limit  of  about  50  eV 
(weakly  dependent  on  plasma  density). 

The  pressure  of  the  collisionless  plasma  in  the  periphery  becomes  anisotropic  as  the 
plasma  contracts  toward  the  axis,  due  to  the  separate  conservation  of  the  electron's  mag¬ 
netic  moment,  and  the  electron’s  angular  momentum  about  the  axis.  Under  these  con¬ 
ditions  the  plasma  can  still  oscillate  in  a  self-similar  manner,  and  for  constant  z-pinch 


current  and  balanced  oscillation  frequencies  (u.’|, , .cl , i*.’*)  the  periphery  can  be  in  an  exact 
equilibrium. 

As  the  current  increases  in  the  power  pulse  part  of  its  energy  is  siphoned  off  by  the 
periphery.  Exactly  what  fraction  is  lost  depends  upon  the  detailed  structure  of  the  low 
density  plasma  corona,  which  is  as  yet  poorly  understood.  For  a  tenuous  periphery,  n 
1013  cm3,  with  the  corresponding  temperature,  T  ~  5.0el  ,  little  power  is  diverted  from 
the  z-pinch.  For  densities  appropriate  to  edge  gas,  n  ~  101S,  cm3,  or  higher  temperatures 
T  ~  50el',the  power  loss  is  substantial. 
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PEOS  Simulations  with  MAGIC. 

Particle  simulations1,2  of  erosion  switches  have  been  notably  divergent  from  the  lab¬ 
oratory  experience  in  several  respects.  The  primary  objectives  of  the  work  described  here 
are:  (i)  to  substantiate,  modify  or  critique  the  earlier  and  contemporary  simulations  (in 
Ref.  1,2)  in  the  context  of  the  DNA  e/m  PIC  code  MAGIC;  (ii)  extend  the  parameter 
domain  of  the  simulations  in  order  to  document  fully  the  implications  of  e/m  PIC  codes  for 
PEOS;  and  (iii)  to  provide  a  comparison  basis  for  evaluating  the  newer  simulation  codes, 
and  novel  algorithms  e.g.  the  magneto-inductive  PIC  code,  or  particle  collision  schemes. 
To  these  ends  a  sequential  upgrade  of  the  PEOS  study  due  to  Grossman1  has  been  under¬ 
taken  with  MAGIC.  The  progression  has  been  (i)  planar  transmission  line  with  H+  ions  at 
1012  cm-3  and  no  radial  drift,  (ii)  cylindrical  geometry,  and  C+  ions  at  1012  cm'3,  and 
(iii)  C+  ions  with  a  radial  drift  ~  2  x  108  cm  sec.  The  planned  continuation  is  to  vary  the 
density  drift-velocity  product,  diagnose  the  axial  momentum  deposition  and  the  electron 
temperature,  and  explore  more  realistic  (BJ5-  like)  downstream  loads. 

A  useful  comparison  of  the  MAGIC  simulation  to  date  with  those  of  Grossman,  et. 
al  (his  fig.  2,3,4)  can  be  drawn  between  the  second  and  third  cases  above.  For  these 
comparisons  the  geometry  of  the  region  is  shown  in  Fig.  5.  The  peak  current  input  to  the 
switch  is  156  k.A  over  5  ns,  and  the  opening  time  (saturated  load  current)  for  case  (ii)  is 
larger  than  10  ns.  This  is  a  somewhat  slower  opening,  and  the  ion  mass  distribution  is 
more  snowplowed  when  the  anode  influx  is  shut  off.  The  current  channel  remains  narrow, 
however,  and  the  anode  end  leads  the  cathode  end.  in  a  manner  similar  to  the  Grossman 
calculation.  A  sequence  of  MAGIC  plots  is  shown  in  Fig.  6,7.  Figure  6  shows  the  rB# 
contours,  while  Figure  7  displays  the  ion  distribution  in  space  at  2.5,  5.0,  7.5,  and  10.0  ns 
during  the  switching  action. 

The  mechanism  of  current  channel  formation  and  switch  opening  shown  here  is  thus 
quite  similar  to  that  depicted  by  Grossman  et.  al  although  the  interpretation  is  subject 
to  some  honest  controversy.  A  central  point  is  the  E  *  B  drift  and  its  relation  to  the  current 
channel,  and  also  the  role  of  ion  drift  in  determining  the  opening  time  and  peak  current.  A 
fair  summary  of  the  case  (ii)  MAGIC  simulation,  in  which  ion  drift  and  regeneration  at  the 
anode  are  absent,  is  that  the  incoming  axial  momentum  of  the  power  pulse  is  coupled  to  a 
small  fraction  of  the  plasma  load  ions  through  a  shock-like  structure  of  charge  separation,  a 
momentum  transfer  front.  At  this  charge  separation,  the  confluence  of  the  charge  density  p 
and  the  E  ■  B  drift  produces  the  current  density  in  the  switch,  c.f.  Fig.  8.  The  opening  is 
forced  when  the  momentum  transfer  front  runs  downstream  of  the  switch  plasma,  causing 
the  active,  charge  carrying  E  ■  B  contours  to  be  only  those  nearest  the  cathode.  One 
can  view  this  process  as  representing  the  completion  of  the  magnetic  insulation  process  in 
the  cathode  sheath.  The  downstream  F,r  still  peaks  near  the  anode  rather  than  near  the 
cathode  as  occurs  upstream,  due  to  the  neutralizing  sheath  space  charge.  The  erosion  of 
ion  density  upstream  of  the  momentum  transfer  front  plays  a  very  weak  role  in  the  process 
here,  when  the  radial  ion  drift  is  set  to  zero:  and  the  opening  time  is  notably  slower  than 
for  the  Grossman  computation. 

When  the  iori  drift)--  2  10*  cm  s)  is  now  turned  on,  case  (iii)  the  opening  time  is 
ml  by  10*  r  to  about  6  ns  as  shown  in  Fig.  0.  10.  This  drifted  1012  cm  3  carbon  plasma  is 
thus  opening  almost  as  last  as  Grossman's  more  slowly  drifting  (1.89  <  108  cm/s)  hydrogen 
plasma  at  a  somewhat  lower  current  level  (lOOkA). 


It  is  common  in  experiment3  to  see  a  wider  current  sheath,  but  the  MAGIC  run  gives 
a  possible  clue  to  a  partial  resolution  of  this  problem  when  current  density  streamline  plots 
and  range  plots  are  examined.  One  finds  circulating  currents  upstream  of  the  momentum 
transfer  front,  such  that  an  axial  B<?  profile  at  fixed  r  would  see  a  slightly  non-monotone 
decay  in  B<?  with  increasing  z.  Experimental  probe  measures  also  hint  at  such  behavior, 
but  no  claim  is  yet  made  as  to  whether  or  not  such  reversals  are  in  fact  resolved  by  the 
experiment.  If  vortices  of  current  density  are  in  fact  present  ,  e.g.  through  gyrokinetic 
modes  upstream  of  the  momentum  transfer  front,  then  an  apparent  broadening  of  the 
current  sheath  would  occur  -  even  though  such  current  is  not  connected  to  the  driving 
circuit. 

The  theory  of  operation  of  PEOS  devices  is  thus  somewhat  confused  and  likely  to 
remain  so  for  a  time.  If  the  e/m  PIC  method  can  show  its  relevance  in  accounting  for 
the  observed  coupling  of  ion  drift  speed  and  opening  time4,  and  if  better  resolution  of  the 
current  channel  structure  can  be  obtained,  then  a  detailed  theoretical  accounting  of  the 
required  collisionality  to  broaden  the  current  channel  properly  would  be  a  logical  step. 
The  means  to  this  assessment  and  the  application  of  more  appropriate  magnetoinductive 
simulation  methods  for  long  time  scale  conduction  phase  operation  are  therefore  among 
our  primary  near  term  goals. 
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I.  Magnetoinductive  Particle-in-Cell  Code  Development 

This  past  year's  report  describes  the  progress  towards  the  utilization  of  Darwin’s1 
magnetoinductive  mode!  for  a  particle  simulation  rode.  The  magnetoinductive  model  is 
desirable  because  it  retains  an  induced  electric  field  while  not  including  light  waves.  T  his 
approximation  relaxes  what  can  be  a  severe  constraint  on  the  time  step  and  must  be  ex¬ 
amined  in  pulsed-power  applications  like  the  VOS.  The  intent  of  any  VOS  is  the  storage 
of  energy  over  a  long  time  frame  (as  measured  in  electromagnetic  signal  transit  times) 
and  then  the  release  of  this  energy  over  much  shorter  times,  more  like  the  electromagnetic 
transit  time,  through  the  interuption  of  a  current  path.  Under  such  conditions  the  dis¬ 
placement  current  cannot  be  neglected  at  the  time  of  opening,  but  it  cannot  be  allowed 
dominance  during  the  conduction  phase  because  it  makes  the  calculation  too  stiff.  The 
Darwin  limit  retains  that  part  of  the  displacement  current  needed  to  resolve  the  interuption 
and  also  removes  the  primary  source  of  stiffness  -  the  traveling  waves. 

However  solving  for  the  magnetic  and  the  electric  fields  becomes  more  difficult  in 
this  limit.  The  magnetoinductive  model  requires  knowledge  of  more  information  at  the 
boundaries.  Substantial  pirogress  has  been  achieved  but  difficulties  still  remain.  The 
next  section  describes  the  magnetoinductive  formulation  and  two  different  schemes  for 
its  implementation.  Results  from  testing  the  model  are  described  in  Section  III.  file 
primary  problem*  encountered  relate  to  the  specification  and  or  convergence  of  boundary 
rondit  ion*. 

II.  The  Magnetoinductive  Formulation 

Part ic|e-in-(  e|l  codes  have  proven  themselves  to  be  valuable  tools  in  examining  prob¬ 
lems  where  the  detail*  of  individual  particle  motion  are  important.  Among  the  most  useful 
and  versatile  of  these  tool*  are  the  fully  electromagnetic  particle  codes  These  codes  fea¬ 
ture  a  complete  classical  description  of  the-  electromagnetics  and  particle  motion.  The 
stabilitv  and  acciiracv  of  the  code*  require  the'  time  step)  to  be  sufficient  Iv  small  to  accu¬ 
ral  el  \  r  O'o|\  e  the  p>t  opaga  t  ion  of  ligli  t  w  ;u  e*.  the  fastest  c  yc  lot  rein  motion,  and  t  he  fastest 
pda-te.a  wave*.  However  'here  a;e  numerous  example*  in  the  pmlsed-power  field  where  '  lie 
vacuum  wavelength*  associate'll  with  a  given  electromagnetic  pulse  are  much  larger  than 
the  maibme  dimension*  hot  these  <a*e*.  Darwin  *  magnetoinductive  model  of  the  elec- 
'  t  a  ami  magnetic  field*  mav  be  u*elul  n.e  Darwin  model  neglect*  retardation  effect*  in 
the  1  ran*  mi-  '  mi  of  mlot  mat  n<n  wfuc  h  tia*  t  e  eflec  t  of  elimmat  ing  purelv  elect  rom  ague!  ie 
mode*  |t  ■*  no  longc  i  eecessarv  to  rest  net  the  time  step  such  that  light  cannot  ptop- 
aga’c  more  •  t  ,  i  r  Mu  •malle’-t  c  c  u  n  pi  i 1  ,•  t  ic  ui  a  I  cell  m  .1  time  tep  \  not  her  advantage  of 

till  •  •  ,  I  [O  I  !  |>  i  •  :  I  I  1 1  1  ■  \  I  o  1 1  I,  1 1  1  I  .  ,1  d  •  .  '  I  I  '  ■  ! 1  1  1 11  i|  \  <  1  I  it  e  |i  -i  I  I  l  >1 1 1  a  g  "  e  •  ■  I  tl  ’  i  I  I  :  I  ,1 1  lot!  .  |  h  e 

i  r  •  1  •  • .  i  i  ,  I  , , .  t  ■  ; ,  1 1  ,  I  ,-*'*•«  i  i  a  . .  |  i.i  *  *  ic  I*  .  1 1 1 1 1  ■  |,|  I  ic  ■  <  i  a,  •  o  ’  ■  ■  c  .i  1 1 1  1 1  n  ,c  b  ■  c  i  I  \  I  a  '  ge 

1 1  i  :  i  e  •  .  •  i  i  a  '  1  i  i ■  i  i  '  ;  ■  •  •  ■  a  g t ; c  ' 1 1  I  d  I  !  1 1  -  c  tl  ■  1 1  1  1 1  a  '  coll  c  a  U  ■  I '  c"  I  M  a  t  e  oV  c  t  '  h' 1 


1  a  '  '  v  ■  i  r  ■ .  •  ii  1 1 


I-  'I.  ■  •  ag-  .  to 


i  Ui  '  o  1 1 1 1  a  g  i  a  a  ic 
■  .  i ;  •  j « f  ■  ■ ' ,  1 1 1 1 , 1 1  1 1 1 1 1 
.1  '  I  t  i  I  (  )  '  I  c  <’t 
:  •  ■  . , .  g 1 1 1  ‘  t c  tic  Id* 


•a  'del .  tie 


ac  Ii  t  line 


'  .  .  [  '  t  <  •  '  •  C  I  I  j  I 


1 1 1  o  I  o  n  g  i '  1 1  c  I  m  a  I 
. .  ■ '  c  I  i  v  •  model 


-.-vv  ys’  --v- 


(2) 

(3) 

(4) 


V  •  E  =  p/(n, 

VxB  =  pnJ  +  \dEi/dt  =  po  Jt< 

V  x  Et  =  ~ dB/dt , 

where  the  transverse  and  longitudinal  fields  by  definition  satisfy 

V  ■  E-p  =  0, 

Vx£l=0. 

This  set  of  equations  automatically  satisfies  the  charge  continuity  equation 

dp/dt  +  V  ■  /  =  0. 

The  field  equations  in  the  magnetoinductive  model  differ  only  in  the  neglect  of  the  partial 
time  derivative  of  the  transverse  electric  field  in  equation  3  from  the  full  set  of  Maxwell’s 
equations. 

The  solution  for  the  transverse  electric  fields  proves  to  be  the  most  difficult  task  for 
the  implementation  of  this  model.  An  explicit  time  differencing  of  Faraday’s  law  leads  to 
violent  numerical  instabilities.  This  time  differencing  tries  to  propagate  information  at  a 
finite  speed  which  is  inconsistent  with  the  instantaneous  Darwin  model.  The  magnetic 
field  is  eliminated  between  equations  3  and  4  to  obtain 


V2Er  —  Mo 


where  for  nonrelativistic  particles  of  species  a: 


(5) 


dJa 

dt 


4 - Ja  x  B  —  <7qV  •  *P a, 
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and  where 


*P a{ x)  =<  nava  va 


>  . 


Inspection  of  equation  6  exhibits  the  difficulty  in  solving  for  the  transverse  electric  field. 
The  time  derivative  of  the  current  is  also  a  function  of  the  transverse  electric  field  that  we 
are  trying  to  solve.  The  combination  of  equations  5  and  6  yields  an  equation  of  the  form. 


V2  Ef  - 


Et  +  Mo 


T 


(7) 


It  is  important  to  notice  that  equation  7  features  a  vector  decompostion  on  the  right  hand 
side  of  the  equation.  This  vector  decomposition  is  analytically  performed  by  a  volume 
integration  of  the  vector  (c.f.  Jackson2).  Equation  7  forms  an  integrodifferential  equation 
that  is  very  difficult  to  solve  directly.  Instead  we  resort  to  iterative  methods. 
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A  standard  way  of  solving  equation  7  suggested  by  Nielson  and  Lewis3  incorporates 
the  subtraction  of  a  constant  times  Ex  from  both  sides.  This  constant  is  chosen  as  some 
mean  value  of  the  term  multiplying  Ex  on  the  right  side  of  equation  7. 

(8) 

T 

where  u>2  is  a  constant.  An  iterative  procedure  is  followed  in  which  the  previous  iteration 

value  of  Ex  is  used  on  the  right  hand  side  of  equation  8  to  find  the  current  iteration  value 
of  Ex-  The  transverse  projection  in  equation  8  is  performed  with  the  aid  of  the  continuity 
equation  in  the  following  way: 
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whence 


dJ_ 

at 
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Boundary  conditions  on  the  second  time  derivative  of  <f>  must  be  specified.  This  proves 
difficult  for  an  arbitrary  configuration.  This  scheme  can  have  problems  converging  if  there 
are  large  variations  in  w2(i)  caused  by  large  density  gradients. 

As  suggested  by  Boyd  and  Hewett4,  an  alternative  formulation  to  equation  8  solves 
the  equation 

V*f  =  «  o^,  ("I 


with 


at 


£  =  Ex  -  Vx- 


Replacement  of  Ex  by  £  +  Vy  in  equation  7  gives 


"?(*) 


w. 


VX  +  PoS, 


(12) 


with 


V2X=  -V-£. 


The  transverse  projection  is  performed  on  the  quantity  ^  to  obtain  Ex  in  this  scheme.  An 
elliptic  equation  is  solved  for  x  to  There  exists  a  great  deal  of  freedom  in  specifying  the 
boundary  conditions  for  Boundary  conditions  on  £  are  determined  by  the  boundary 
conditions  on  Ex  and  on  the  first  and  second  derivatives  of  x  normal  to  the  boundary.  A 
disadvantage  is  that  the  boundary  conditions  om  $  are  changed  with  each  iteration. 


III.  Tests 
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The  magnetoinductive  code  that  has  been  developed  is  benchmarked  in  a  series  of  tests 
described  in  this  section.  The  first  such  test  which  directly  compares  the  two  schemes  is 
presented  first.  An  initially  quasineutral  plasma  with  an  equal  number  of  protons  and 
electrons  is  placed  in  a  box.  The  box  is  positioned  in  an  azimuthally  symmetric  cylin¬ 
drical  geometry.  The  walls  of  the  box  are  assumed  to  be  conducting  and  are  set  to  zero 
potential.  Particles  within  the  box  elastically  rebound  off  the  walls.  There  are  not  any 
external  sources  of  energy  so  there  should  not  be  any  energy  flow  in  or  out  of  box.  The 
plasma  is  initialized  in  an  equilibrium  configuration  where  the  plasma  density  is  uniform 
and  the  electron  and  ion  temperatures  are  equal  and  constant.  For  the  test  presented 
here  the  plasma  is  placed  in  the  center  half  of  the  box  with  vacuum  in  the  remainder  as 
approximately  shown  in  Fig.  1.  The  inner  box  radius  is  2.5  cm  and  and  the  outer  radius 
is  5.0  cm.  The  box  is  10.0  cm  long  in  the  z-direction.  The  grid  is  20  cells  by  20  cells.  The 
plasma  is  assigned  a  temperature  of  1  keV  and  the  plasma  density  is  1.0  x  1012  cm-3. 
The  time  step  is  4.0  x  10~12  seconds. 

The  results  from  Nielson  and  Lewis  scheme  are  physically  reasonable.  The  simulation 
is  followed  forward  in  time  for  2.0  ns.  Electrons  expand  faster  into  the  vacuum  region 
because  their  thermal  velocities  are  larger.  Local  charge  imbalances  are  created  with  a 
negative  charge  density  occurring  in  what  was  originally  the  vacuum  and  a  positive  charge 
density  occurring  in  the  main  body  of  the  plasma.  The  contours  of  the  electrostatic 
potential  in  Fig.  2  highlight  this  behavior.  The  maximum  of  the  electric  potential  is  in 
the  same  range  as  the  plasma  temperature.  It  is  unlikely  that  potentials  much  larger  than 
the  temperature  would  become  established  since  the  plasma  is  initialized  in  a  thermal 
equilibrium.  Towards  the  beginning  as  electrons  stream  out  of  the  plasma,  a  current 
towards  both  walls  of  constant  z  builds  up.  This  creates  a  buildup  of  magnetic  flux.  Later 
in  time,  electrons  that  hit  the  wall  are  reflected  reversing  their  direction  of  motion.  As 
these  electrons  come  counterstreaming  back,  the  net  current  starts  decreasing  and  the 
magnetic  flux  decreases.  This  decrease  in  the  magnetic  flux  is  seen  in  the  plot  Err  in  Fig. 
2.  Considering  the  coarseness  of  the  grid,  it  is  seen  in  Fig.  3  that  energy  is  quite  well 
conserved  increasing  by  only  0.7%  over  the  length  of  the  run.  In  runs  where  the  plasma 
completely  filled  the  computational  volume  and  the  other  parameters  remained  constant, 
energy  is  conserved  to  0.2%. 

In  contrast  the  formulation  of  Boyd  and  Hewett  was  unstable  quickly  and  stopped 
after  a  few  tens  of  time  steps.  The  results  in  Fig's.  4  and  5  show  informations  after  10 
time  steps  or  .04  ns.  In  Fig.  4,  the  plots  of  the  electrostatic  potential  have  the  wrong 
sign  in  the  vacuum  and  plasma  regions.  The  maximum  value  of  30  kilovolts  is  much 
larger  than  what  is  expected.  Similarly  the  maximum  value  of  the  transverse  electric 
field  is  much  larger.  As  can  be  seen  from  Fig.  5  the  conservation  of  energy  is  badly 
degraded  in  this  example  after  the  simulation  has  progressed  forwards  only  10  time  steps. 
This  scheme  is  clearly  unavceptable  for  this  particular  problem.  It  should  be  mentioned 
that  this  is  probably  the  worst  type  of  problem  for  this  scheme.  Its  worst  failing  is  that 
the  boundary  conditions  for  the  auxiliary  function  \  are  changed  with  each  iteration  of 
the  solution.  In  this  case  the  fields  have  random  character  that  are  determined  by  the 
motion  of  the  particles.  The  fluctuations  leave  ET  having  a  spiky  and  random  character 
that  can  change  considerably  from  time  step  to  time  step.  The  iterative  procedure  has  a 
great  deal  of  difficulty  in  converging  the  boundary  conditions  in  those  circumstances.  This 
technique  is  probably  suitable  in  a  situation  where  external  conditions  primarily  determine 


2.  Contours  of  the  electrostatic  potential  and  transverse  electric  field  at  2  ns 
Nielson  and  Lewis  scheme  for  the  configuration  shown  in  Fig.  1.  Maximum 
ire  \<t>\maz  ~  1.5  kV  and  \ET\max  =  1  kV/M. 
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nservation  for  the  Nielson  and  Lewis  scheme  for  the  configuration 
■  evolving  for  2  ns  Maximum  change  is  0.7%. 


4.  Contours  of  the  electrostatic  potential  and  transverse  electric  field  at  2  ns 
Nielson  and  Lewis  scheme  for  the  configuration  shown  in  Fig.  1.  Maximum 
ire  I <f>\max  =  30  kV  and  |J5rjmoi  =  70  MV/M. 
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are  imposed.  The  solution  for  the  electrostatic  potential  is  given  by 
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The  current  can  be  determined  by  eliminating  the  magnetic  field  from  Eq’s  15a  and  16  to 
get, 

d[a(z ,0  2  n  (q  d<f> 

dz  fn(£)  dt  ' 

which  has  the  solution 

T  (n  y  t\  -  2  ^  €°  ( Zl  ~  Z)  ^  T  <M\ 

4  ’  ’ 0  "  ln(i)  dt  h'  (17) 

where  1 1  is  the  current  flowing  through  the  resistor.  Thereby  Be  is  given  by 


Be  =  - 


2  7T  r 


(Zl  -  z)  d<fr 
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Lastly  the  solution  for  the  transverse  part  of  the  electric  field  is  given  by  the  z- 
component  of  Faraday’s  Law  with  equation  18: 


dETr 

dz 


Ho  dIL  (zi  -  z)  d2<f> 

2  n  r  dt  r  c2  /n(J)  a<2  ‘ 


Since  the  load  is  assumed  to  be  purely  resistive,  we  have  as  a  boundary  condition  at  the 
load  side 

ETr{r,  zi,t)  =  0. 

Integrate  once  with  respect  to  z  with  the  above  boundary  condition  to  get 


ETr(r,z,t)  =  - 
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The  agreement  between  the  analytic  solutions  given  by  Eq’s  18  and  19  and  the  code  results 
is  excellent  as  can  be  seen  in  Fig.  7.  The  results  in  Fig.  7  were  obtained  with  a  20<  20 
grid.  This  test  demonstrates  that  the  routines  that  solve  for  the  electric  and  magnetic 
fields  are  working  quite  well. 

1  he  next  test  concentrates  on  the  particle  motion  and  emission.  The  results  from  the 
simulation  are  compared  with  the  one-dimensional  C'hild-Langmuir  diode  in  an  electro¬ 
static.  nonrelativist ic  limit.  The  geometry  for  this  test  is  shown  in  Fig.  8.  The  magnetic 
field  and  the  inductive  electric  held  are  removed  from  the  simulation.  The  only  field  that 
effects  particle  motion  is  the  electrostatic  electric  field.  A  positive  potential  is  placed  on  the 
top  plate  with  respect  to  the  bottom  plate.  Electrons  are  emitted  from  the  bottom  plate 
arid  travel  to  the  top  plate  where  they  are  absorbed.  It  should  be  mentioned  that  for  this 
test  electron  emission  is  radial  ami  the  equations  presented  here  for  the  Child-Langmuir 
diode  m  a  <  artesian  geomet  ry.  However  if  the*  hange  in  radius  from  t  he  bottom  to  t  he  top 
is  small  <  orupared  to  the  average  radius,  the  emission  should  be  approximately  Cartesian. 
Fieri  mu  How  from  the  bottom  pint**  is  presumed  to  be  space-charge  limited. 

I  he  <  hild- Langmuir  diode  finds  t  he  stationar\  solutions  when  electrons  are  assumed 
to  he  <  old.  i  e  fin  it  e  t.emperat  ure  effer  t  s  are  ignored .  The  voltage  V  ( t )  of  the  top  plate  wit  h 


Figure  7.  Compares  the  analytic  solutions  with  tjie  code  solutions  in  a  vacuum. 
Agreement  is  very  good. 
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Such  a  splitting  depends  on  t  he  external  circuits  that  the  simulation  is  connected  to.  This 
area  is  under  investigation  at  this  time.  Another  possibility  is  to  artificially  limit  in  a  wav 
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Figure  9.  Simulation  results  favorably  reproduces  potential  expected  for  Child 
Langmuir  diode. 


that  dampens  the  instability.  This  approach  may  be  useful  in  achieving  a  time  stationary 
solution  but  would  not  be  realistic  enough  to  follow  the  time  evolution  of  a  problem. 
The  magnetostatic  code  ELECTRA  is  probably  more  suitable  for  these  applications.  The 
next  steps  are  to  include  more  realistic  circuit  models.  These  would  be  added  first  to  an 
electrostatic  code  and  subsequently  implemented  in  magnetostatic  and  magnetoinductive 
codes. 

IV.  Summary  and  Expected  Applications 

The  magnetoinductive  formulation  is  promising  in  the  amount  of  physics  that  is  in¬ 
cluded  while  not  being  restricted  to  resolving  propagation  of  light  waves.  Practical  im¬ 
plementation  has  proved  to  be  difficult.  All  of  the  various  pieces  of  the  code  have  passed 
benchmark  tests.  The  primary  difficulty  lies  in  trying  to  specify  the  boundary  condi¬ 
tions.  Two  different  schemes  have  been  examined  for  performing  the  iterative  solution 
on  the  transverse  electric  field.  The  Nielson  and  Lewis  scheme  is  preferred  but  requires 
specification  of  the  potential  associated  with  tee  longitudinal  electric  field  and  its  first- 
and  second-order  time  derivatives.  This  is  not  an  easy  task  in  an  arbitrary  configuration. 
The  Boyd  and  Hewett  scheme  requires  less  information  on  the  boundaries  but  has  more 
problems  converging  when  the  transverse  electric  field  is  not  sufficiently  smooth  at  the 
boundaries.  The  utility  of  the  magneoinductive  formulation  depends  on  the  particular 
application. 

The  initial  applications  of  the  magnetoinductive  codes  are  aimed  towards  the  fast 
opening  switches  crucial  to  inductive  energy  storage  schemes.  These  will  emphasize  looking 
at  important  pieces  of  the  problem  in  a  simpler  configuration.  One  such  application  is  to 
examine  the  role  microinstabilities  play  in  the  current  conduction  phase  of  PEOS.  Another 
is  to  implement  different  models  of  particle  emission  and  investigate  how  these  different 
models  effect  current  conduction  throught  the  sheath.  It  may  also  be  used  to  examine  the 
reflex  switch. 
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Numerical  Models  of  Bremsstrahlung  Diodes 

Much  effort  in  collaboration  with  Rerkeley  Research  Associates  during  the  past  year 
has  been  aimed  at  improvements  in  modeling  Bremsstrahlung  diodes.  Bremsstrahlung 
diodes  are  the  most  common  x-ray  source  in  use  today.  Although  their  spectrum  is  some¬ 
what  harder  than  that  desired  b\  many  users,  their  simplicity  and  reliability  as  a  radiation 
source  has  made  them  popular.  The  trend  has  been  for  x-ray  source  users  t  >  demand  longer 
pulses  and  higher  overall  dose  over  a  broader  area  while  maintaining  acceptable  uniformity. 
This  trend  has  pushed  present  diode  designs  to  their  limits  leaving  their  performance  in 
this  regime  somewhat  uncertain. 

The  technical  issues  that  arise  in  this  new  diode  performance  regime  are  similar  to 
those  that  have  arisen  in  the  past.  In  particular,  diode  design  requires  a  detailed  picture  of 
electron  flow  from  the  cathode  to  the  converter  as  well  as  global  electrical  characteristics 
such  as  impedance.  In  addition  to  these  issues,  the  longer  pulse  times  bring  other  effects 
into  consideration  such  as  gap  closure  from  plasma  formation  at  the  cathode  and  anode 
walls. 

In  cases  where  theoretical  support  of  diode  design  has  been  helpful,  t  lie  results  have 
come  mainly  from  electromagnetic  particle-  in-cell  (PIC,')  simulations.  Such  simulations 
follow  the  trajectories  of  representative  electrons  while  advancing  the  electromagnetic  field 
equations  on  a  grid  of  cells  in  a  way  that  is  consistent  with  the  currents  created  by  the 
electrons  themselves.  Unfortunately,  such  techniques  have  also  encountered  limitations  as 
pulse  times  have  become  longer.  Technically  to  advance  the  electromagnetic  field  stably, 
the  timestep  must  be  substantially  smaller  than  the  time  it  takes  for  a  light  wave  to  cross 
a  cell  in  the  calculation.  In  low-  impedance  diodes  where  the  gap  spacing  is  a  centimeter 
or  less  and  at  least  ten  or  twenty  cells  are  needed  to  span  the  gap  this  would  require 
sub-picosecond  timesteps.  With  pulse  times  of  hundreds  of  nanoseconds,  the  number  of 
timesteps  for  a  complete  simulation  could  run  into  the  millions.  This  kind  of  numeri¬ 
cal  difficulty  in  which  timescales  that  must  be  resolved  are  very  far  apart  is  known  as 
"stiffness" . 

Among  the  aims  of  our  current  research  is  to  find  useful  numerical  techniques  ior 
diode  modeling  that  will  work  for  these  longer  pulse-lengths.  In  similar  circumstances 
where  very  fast  timescales  interfere  with  the  slower  timescales  of  interest,  two  general 
strategies  are  used.  The  first  is  to  try  a  reduced  physical  description,  viz.  a  simpler  form 
of  the  governing  equations  in  which  the  fast  timescale  is  absent.  The  magnetoinduct  ive 
model  and  the  steady  state  model,  Plectra,  are  examples  of  this  approach.  The  second 
strategy  is  to  find  other  numerical  formulations  of  the  complete  equations  in  which  'ast 
oscillations  are  stably  damped  away  when  the  timesteps  are  too  large  to  resolve  them. 
This  approach  is  being  investigated  as  well. 

Steady  state  diode  calculations  and  development  of  ELECTRA 

The  long  pulse  lengths  of  interest  in  many  applications  often  mean  that  the  power  flow 
in  the  diode  lias  readied  a  steady  state.  Keen  in  circumstances  where  this  is  not  precisely 
true,  the  steady  state  or  equilibrium  picture  of  the  diode  represents  the  main  trend  in  its 
behavior  and  provides  insight  that  can  be  exploited  in  future  designs.  The  electromagnetic 
field  equations  appropriate  to  stationary  charge  distributions  and  currents  are  simply: 


Poisson's  Equation,  and 

V  1  B  /inj. 

Ampere’s  Law,  with 

E  -  Vo  and  V  B  0. 

Here  E  is  the  electric  field,  and  B  is  the  magnetic  induction,  o  is  the  static  potential. 
J  is  the  current  density  and  p  is  the  charge  density.  These  field  equations  can  he  used  as 
the  basis  of  a  particle-in-cell  simulation  in  which  the  particles  are  moved  a  timestep  and 
then  the  charge  and  current  distributioi  s  are  formed  from  them  and  taken  as  statu  (or 
the  purpose  of  calculating  the  fields.  Assuming  a  sufficient  number  of  timesteps  and  t  he 
existence  of  a  steady  state,  the  currents  and  charge  distributions  indeed  become  steads 
enough  for  the  static  field  assumption  to  be  valid. 

The  static  field  assumption  offers  two  advantages  in  the  present  application.  Since 
electromagnetic  waves  do  not  exist  in  t he  st at  ic  limit .  there  is  no  t  imestep  const  rain t  aris  ing 
from  them.  While  accuracy  demands  that  the  particles  themselves  not  move  more  than  a 
cell  per  t  imestep,  this  is  a  milder  limit  at  ion.  Secondly,  since  the  steady  state  sol  ut  ion  is  the 
aim  of  the  calculation  and  because  the  specific  intermediate  state's  are  not  critical  to  the 
attainment  of  a  solution,  any  expedient  measure  taken  to  speed  convergence  to  equilibrium 
is  allowed.  In  short,  it  should  be  possible  to  arrive'  at  a  steadv  slate-  solution  in  fat  fewer 
timesteps  using  the  static  equations  than  by  using  the  full  set  of  Maxwell's  equations 

Simulations  based  on  static  field  equations  have  been  use-d  in  the-  past  at  the  Naval 
Research  Laboratory.  Sandia  National  Laboratory  and  elsewhere.  They  represent  some  of 
the  earliest  attempts  at  diode  mode-ling:  e  g.  the  DIODE20  code  at  NRL.  They  arc  not 
common  today  for  two  reasons:  the  high-impedance,  short  pulse  diodes  ran  fie  mode  led 
with  the  more  complete  elect  romaget  ic  simulations;  and  t lie-re-  is  no  general  purpose  easilv 
used  static  code  available  to  the  pulsed  powe-r  community  as  is  the  case  with  the  fu!l\ 
electromagnetic  codes  (e.g.  MAOIC.  MASK). 

Berkeley  Research  Associates  in  collaboration  with  NRL  has  revived  tin-  equilibrium 
diode  code  for  long  pulse  diode  applications.  KLKCTRA  is  the  name  of  the  new  code 
produced  under  this  effort  and  test  versions  have  been  in  application  over  the  past  six 
months.  KLKCTRA  has  features  in  common  with  the  old  DIODK2D  code  once  usee  at 
NRL.  It  is  a  two-dimensional  simulation  in  r  an  /  allowing  anode  and  cathode  shape-  to 
be  represented  as  internal  structures  conformal  to  the  grid 

Simple  field  emission  algorithms  are  based  on  Causs'  Law.  Both  KLKCTRA  and 
DIOPK2I)  have  simple  foil  converter  models  for  reflex  diodes.  KLKCTRA.  however,  lias 
many  improvements  over  the  older  code. 

KLKCTRA  is  based  on  optimized  fini't  difference  field-solvers  rather  than  Fourier 
transform  methods.  Particle  motion  is  referred  to  a  "grid-based"  coordinate1  system  rather 
than  actual  physical  coordinate's.  This  combi  nat  ion  permits  non  uniform  gt  id  ding  in  eit  her 
coordinate  with  arbitrary  spacing. 

The  magnetic  field  calculation  in  KLKCTRA  can  be  based  either  on  instantaneous  or 
time-averaged  current  distributions  which  can  speed  convergence  to  steadv  state. 

Improved  particle  movers  are  incorporated  into  KLKCTRA.  One  option  allows 
timesteps  much  larger  than  the  cyclotron  period  (HA t  1).  The  methode  uses  a  magnetic 
moment  correction  which  is  inoperative  when  12 A /  v  1,  but  produces  the  correct  E  B 
drifts  when  f2Af  :>»  1.  This  improvement  allows  simulation  of  pinch-beam  diodes  without 
timestep  limitations  imposed  by  the  gyrofrequency. 
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Suhc\  cling,  of  fast  relat  ivist  k  cle<  I  nuts  tela'  i\ e  to  t  he  slower  elect rons,  and  of  elect  rons 
vs.  ions  will  he  incorporated  to  speed  calculation. 

We  will  experiment  with  dynamic  regridding  or  mult i-gridding  to  reduce  the  overall 
number  of  timesteps  m  the  calculations. 

Some  experiments  with  curvilinear  coordinates  and  finite  element  field  calculations 
will  he  undertaken.  These  ma\  lead  to  improvements  in  the  realism  of  cathode  and  anode 
shapes  and  also  give  better  field  emission. 

Applications  of  ELECTRA 

Although  primarily  in  'he  development  stage  in  ItWi,  El.FiCTRA  has  hern  apolied  to 
diode  design  in  one  instance  and  is  now  finding  application  in  several  others. 


Inverst*  Reflex  Diode  application. 

I  he  ti rs t  application  of  the  prototype  ELECTRA  early  this  year  was  to  the  Inverse 
Wefley  Diode  1  Figure  1  below  illust rates  a  typical  design  with  the  electron  flow  channeled 
along  the  axis  and  then  directed  outward  in  radius.  This  is  distinct  from  the  pinch  beam 
diodes  in  which  the  power  density  increases  with  decreasing  radius.  The  inverse  diodes 
nave  been  suggested  as  a  benign  load  for  pulsed  power  experiments  as  well  as  an  advanced 
(onwrter  design  for  broad  area  illumination.  Simulations  with  ELECTRA  confirmed  the 
tendeiM  \  of  ion  emission  to  complicate  and  sometimes  spoil  the  diode’s  performance  as 
was  suspected  after  preliminary  experiments.  ELECTRA  provided  the  asymptotic  state 
of  the  diode  giving  insight  into  its  dynamical  behavior.  Because  of  its  flexibility  over  the 
older  static  codes,  several  designs  could  be  simulated  in  a  short  time. 

Blackjack  ring  diode  simulations 

The  new  demands  of  pulsed  power  include  higher  dose  as  well  as  better  dose  unifor¬ 
mity  One  successfnl  diode  design  in  this  regard  has  been  the  modular  ring  diodes  on  the 
Blackjack  simulator  at  Maxwell  Laboratories.  A  schematic  representation  of  the  diode  as 
seen  from  the  back  is  shown  in  figure  2.  In  this  design  the  electric  pulse  from  the  generator 
follows  the  posts  located  around  the  ring  through  a  hole  in  the  anode  to  the  azimuthally 
symmetric  ring  cathode.  The  cathode  is  notched  as  shown  in  cross  section  (figure  3)  pro¬ 
ducing  two  electron  beams.  These  beams  impinge  on  the  anode  and  the  return  paths 
divide  according  to  the  inductance  of  the  transmission  line  producing  current  above  and 
below  the  beams.  The  magnetic  field  pressure  produced  on  the  sides  of  the  beams  forces 
pinching  toward  the  center  of  the  anode. 

The  virtue  of  modular  ring  diodes  lies  in  the  ability  to  place  rings  concentrically  wit  bin 
ot her  rings  to  form  arrays.  Such  an  array  offers  higher  dose  than  a  single  converter.  Several 
technical  issues  remain  before  this  design  can  be  exploited  to  the  fullest  advantage,  and 
simple  questions  arise  for  which  simulations  may  give  answers.  For  example,  what  an  *b« 
different  es  between  the  electron  flow  in  rings  at  small  radii  versus  those  at  larger  ram 
Or  put  another  way,  how  critical  is  the  ring  geometry  as  a  function  of  radius'.’  Fm  la'g* 
arrays  to  work  well  each  ring  must  perform  properly.  The  basic  electrical  properties  <d  •  f  • 
rings  (e.g.  impedance)  must  also  be  well  understood  so  that  circuil  models  for  tin  ,rg< 
diode  arrays  will  be  accurate.  At  present  the  closest  packing  of  ring  diodes  in  radi  o 
about  10cm.  Designers  would  like  to  pack  diodes  closer  to  provide  latter  unitor  m  ••  •.  : 

one  needs  to  know  how  close  the  anode  can  come  to  the  cathode  shank  befon  -•  df 
insulation  breaks  down.  Another  question:  how  thin  can  the  cathode  be  mad* 
close  can  the  beams  be  brought  together)  without  affecting  performance’ 

We  have  begun  to  use  ELECTRA  to  answer  some  of  these  que«t ;<.!,<  \ 

preliminary  run  is  shown  below  in  figure  d.  Here  we  see  electron- cm  •• 

ring  cathode  and  pinching  at  the  anode.  Such  plausible  te-ult-  r*  <  •  ■  k  , 
encouraging  is  the  small  amount  of  computer  time  ( ies-  tim  • 

783  at  NR L)  needed  to  obtain  the  solution.  Berkeley  Re-.  ■  \ 

this  study  in  collaboration  with  Branch  1720  <:  \Ri  .,n 
Laboratories.  We  hope  to  be  Mde  to  pure  '-.get',,  -  >...*> 
results  to  help  improve  t  his  ronvi  :  ter  u  •  <  . 
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The  radial  implosion  of  hollow,  cylindrical  neon  gas  columns,  driven  by  currents  of  up  to  1.45 
M  A,  produces  a  linear  Z  pinch  with  over  70%  of  the  radiation  in  neon  K  lines.  A  plasma  erosion 
opening  switch  (PEGS)  is  used  to  eliminate  prepulse  ar.d  to  reduce  the  current  rise  time  from 
—60  to  ~20ns.  Incorporation  of  the  PEGS  improves  the  uniformity  of  the  Z  pinch  and  increases 
the  radiation  yield. 


Pulsed  power  generators  have  been  used  since  the  early 
1970's  to  implode  cylindrical  plasma  columns  for  the  pro¬ 
duction  of  hot,  dense  Z  pinches.  This  research  has  evolved 
from  the  explosion/implosion  of  a  single  wire, 1  to  the  implo¬ 
sion  of  cylindrical  arrays  of  wires2  and  annular  gas  jets.3  In 
this  letter,  we  present  the  first  published  experimental  re¬ 
sults  on  the  implosion  of  annular  neon  gas  jets,  both  with  and 
without  the  incorporation  of  a  plasma  erosion  opening 
switch  (PEOS)4  to  modify  the  driving  current  shape. 

The  Naval  Research  Laboratory’s  Gamble  II  generator 
can  deliver  up  to  2  TW  of  electrical  power  to  a  2-fl  load  in  60 
ns.  In  the  experiments  described  here,  up  to  1.45  MA  was 
delivered  to  a  —  50  nH  gas-puff  load  prior  to  implosion  with 
no  PEOS  [Fig.  1(a) ),  Slightly  over  1  MA  was  delivered  to 
85-nH  load  when  the  PEOS  was  used  [Fig.  1(b)  J.  An  annu¬ 
lar  gas  puff,  2.5  cm  in  mean  diameter,  was  formed  with  a  gas 
injection  system  and  supersonic  nozzle  provided  by  Physics 
International  Company.5  The  neon  gas  was  preionized  with 
a  UV  dashboard  source  a  few  (is  before  the  current  pulse.  A 
thin  wire  mesh  anode,  4  cm  from  the  nozzle,  allowed  the  gas 
to  flow  into  a  dump  tube.  The  current  returned  via  metallic 
posts  at  a  diameter  of  7  cm  allowing  90"  diagnostic  viewing. 
Figure  1(a)  shows  the  current  with  no  PEOS  and  the  result¬ 
ing  neon  X -shell  radiation  pulse  obtained  with  a  filtered  x- 
ray  diode.  The  delay  of  the  x  rays  relative  to  the  current 
onset  is  1 10  ns.  The  equivalent  current  and  x-ray  traces  with 
a  PEOS  are  shown  in  Fig.  1(b).  The  rise  time  of  this  current 
is  only  20  ns.  The  delay  of  the  x  rays  relative  to  the  current  is 
the  same  as  in  Fig.  1  Ca).  Thus,  the  time  required  to  implode 
the  plasma  is  approximately  the  same  with  and  without  the 
PEOS 

The  PEOS  technology  utilized  is  a  recent  development 
and  is  the  subject  of  ongoing  research.*  Suffice  to  say  that  the 
PEOS  initially  serves  as  a  short  circuit  to  the  generator  pulse. 
When  the  switch  current  reaches  a  predetermined  value,  it 
rapidly  opens  diverting  the  current  to  the  load. 

Absolute  x-ray  yields  were  measured  with  filtered  x-ray 
diodes,  tantalum-foil  calorimeters,  bolometers, *  and  curved 
crystal  spectrographs  Total  yields  of  neon  X-shell  radiation 
(near  1  keV)  measured  with  these  diagnostics  were  consis¬ 
tent  among  themselves  to  within  £  15%.  These  yields  were 
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measured  as  a  function  of  plenum  pressure  for  five  driving 
current  levels  with  no  PEOS  [Fig.  2(a)].  The  curves  are 
upper  bounds  of  the  data  scatter  for  each  current  indicated. 
Typically  a  -f  0%,  —  20%  uncertainty  applies  to  each 
curve  due  to  shot-to-shot  variations  in  the  peak  current  and 
for  gas-puff  variations  due  to  nozzle  degradation  in  repeated 
shots.  For  each  driving  current,  the  yield  peaked  at  a  mass 
per  unit  length  ( m )  estimated  by  thin  shell  calculations1  to 
vary  between  —  1 5  and  45  /rg/cm  at  the  lowest  and  highest 
currents,  respectively.  On  some  shots,  an  unfiltered  bolom¬ 
eter6  was  employed  simultaneously  with  the  filtered  diag¬ 
nostics.  These  measurements  indicate  a  production  effi¬ 
ciency  of  A -shell  radiation  which  exceeds  70%  of  the  total 
radiated  energy.  Figure  2(b)  shows  a  plot  of  the  five  peak 
yields  from  Fig.  2(a)  as  a  function  of  mass  and  peak  current 
( I) .  The  strong  scaling  (Ka/Vm1)  indicated,  is  consistent 
with  1-D  slug  model  calculations7  and  results  from  other 
laboratories." 

X-ray  yield  measurements  versus  pressure  are  more  dif¬ 
ficult  with  the  PEOS  since  both  the  shape  and  magnitude  of 
the  current  depend  critically  on  the  injected  PEOS  plasma 
flux  and  on  the  timing  of  plasma  injection  relative  to  current 
arrival.  Shot-to-shot  variations  in  these  quantities  lead  to 
load  currents  of  different  shapes  and  amplitudes,  thus  inhi¬ 
biting  the  systematic  variation  of  gas-puff  pressure  at  a  near¬ 
ly  constant  load  current.  This  problem  is  particularly  acute 
due  to  the  strong  dependence  of  AT-shell  radiation  on  current 
and  mass.  Nevertheless,  X-shell  yields  have  been  measured 
with  the  PEOS  for  a  range  of  current  values  and  the  results 
are  shown  in  Fig  3.  Each  solid  dot  represents  a  shot  with  the 
PEOS  ignoring  the  mass  or  the  exact  shape  of  the  current 
For  all  these  shots,  the  current  rose  in  10-30  ns,  significantly 
less  than  the  60-ns  rise  time  without  the  PEOS.  These  data 
are  compared  to  the  optimum  results  without  the  PEOS, 
represented  by  Ihc  solid  curve.  In  general,  yields  with  the 
PEOS  are  higher  than  the  optimum  yields  without  the  PEOS 
at  a  given  current.  Some  shots  show  an  increase  in  yield  of  as 
much  as  a  factor  of  2  5  (dashed  line) 

Images  of  total  X-shell  emission,  obtained  with  time- 
integrated  pinhole  photography  and  with  spatially  resolved 
x-ray  spectroscopy,'*  show  qualitative  differences  with  ver¬ 
sus  without  the  PEOS  ( Fig.  4 ) .  They  show  the  Z  pinch  to  be 
more  uniform  axially  and  tighter  radially  when  the  PEOS  is 
used  For  example,  the  image  is  partially  obscured  by  one  of 
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FIG  I  Arrangement  of  the  gas  puff  on  Gambit  II  (a)  without  the  PEOS  and  (b)  with  the  PEOS.  Load  current  and  K  shell  x-ray  traces  are  shown  for  each 

case. 


PLENUM  PRESSURE  (psiA) 


PEAK  LOAD  CURRENT  IMA) 


the  return-current  rods  without  the  PEOS,  while  it  is  fully 
visible  with  the  PEOS.  Although  the  sausagelike  character 
of  the  pinch  is  still  evident  in  the  latter,  plasma  flaring  is  less 
severe  and  is  superimposed  on  an  intense  core  <  4  mm  in 
diameter  which  is  reasonably  constant  along  the  z  axis.  The 
axial  uniformity  of  the  pinch  cannot  be  easily  judged  in  these 
photographs  due  to  overexposure.  However,  spectral  lines, 
spatially  resolved  along  the  pinch  axis,  clearly  show  a  de¬ 
crease  in  intensity  from  the  cathode  to  the  anode  in  the  non- 
PEOS  case.  In  addition,  the  line  intensity  ratio  varies  notice¬ 
ably  along  the  axis.  The  Ly-a/He-a  ratio  decreases  from 
~0.9  at  the  cathode  to  —0.2,  2  cm  from  the  cathode.  In 
contrast,  for  the  PEOS  case,  the  line  intensities  are  nearly 
constant  along  the  pinch  axis  and  the  ratio,  Ly-a/He-a,  is 
2  1.1  along  the  length  of  the  pinch.  A  collisional  radiative 
equilibrium  (CRE)  model10  was  used  to  estimate  the  plas¬ 
ma  parameters  consistently  with  the  experimental  results. 
The  model  matched  the  average  Ly-a/He-a  ratios,  radiation 
times,  and  total  yields  observed  experimentally  but  failed  to 
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FIG  ?  Neon  A -shell  yields  measured  withoul  the  PLOS  <  a  )  Yield  vs  gas  FIG  1  Neon  A -shell  yields  with  the  PEGS  as  a  function  nf  peak  load  cur 

pulT  plenum  pressure  for  different  peak  load  currents  and  (b)  peak  yield  rent  The  solid  curve  is  Ihc  optimum  yield  without  (he  PEGS  scaled  as  /‘ 

scaling  with  [>cak  load  current  and  plenum  pressure  The  dashed  curve  is  2  5  times  the  solid  curve 
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FIO.  4.  iC-ahell  x-r*y  pinhole  photographs  (top)  and  spatially  resolved  Af- 
hne spectra  (bottom)  for  neon  unpkajonj  without  and  with  the  PEOS.  The 
cathode  nozzle  is  on  the  left  and  the  wire  mesh  anode  is  on  the  nght. 


match  the  higher  order  line  ratios  ( Ly-a/l.y-0,  H.e-a/Hc-0, 
etc.).  This  inconsistency  is  not  presently  understood  and  is 
under  investigation.  The  Ly-a  and  He-a  lines  are  dominant, 
with  ~90%  of  the  total  radiated  line  energy,  making  this 
plasma  a  nearly  monochromatic  intense  x-ray  source  near  I 
keV.  Based  on  the  axially  averaged  Ly-a/He-a  ratios  we 
estimate  a  plasma  temperature  of  150  eV  without  the  PEOS 
and  185  eV  with  the  PEOS.  These  numbers  are  consistent 
with  loadings  of  34  and  29  /rg/cm,  respectively,  and  a  1.8- 
mm-diam  radiating  core  for  both.  We  conclude  that  both  a 
hotter  and  longer  radiating  plasma  contribute  to  the  ob¬ 
served  increase  in  AT-shell  radiation  with  the  PEOS. 

In  summary,  neon  gas  puffs  driven  with  1.25  MA  from 
the  NRL  Gamble  II  generator  emit  more  than  70%  of  the 
total  radiated  energy  in  the  K  shell.  This  yield  scales  as  the 
fourth  power  of  the  current  with  ~4  kJ  being  radiated  at 
optimized  mass  and  1  45  MA  of  current.  For  a  given  peak 
current,  the  K  shell  yield  can  be  increased  by  as  much  as  2.5 
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by  reducing  the  current  rise  time  with  a  PEOS.  Use  of  a 
PEOS  also  improves  pinch  quality  in  the  radial  direction  and 
its  axial  uniformity.  We  speculate  that  the  elimination  of  the 
early  part  of  the  current  with  the  PEOS,  limits  the  growth 
time  of  instabilities  resulting  in  a  more  orderly  collapse  of  the 
cylindrical  plasma  on  axis. 
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